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ABSTRACT 


The  m i c rop i pe t te  capillary  electrode  is  being  used 
extensively  in  current  physiological  research  in  the  meas¬ 
urement  of  the  action  potentials  of  nerve  and  muscle  cells. 
Electrically,  the  electrode  is  equivalent  to  a  distributed 
low-pass  R-C  transmission  line  which  will  distort,  in  a 
definite  way,  any  applied  potential.  In  the  past,  based  on 
a  very  approximate  equivalent  circuit,  the  generation  of 
negative  capacity  at  the  amplifier  input  has  been  used  to 
approximately  compensate  for  this  distortion.  However,  the 
best  negative  capacity  adjustment  differed  from  an  exact 
neutralization  of  the  input  capacity.  This  investigation 
is  concerned  with  a  compensation  scheme  employing  the  prin¬ 
ciples  of  network  equalization  and  a  distributed  R-C  network 
which  will  allow  a  signal  entering  the  m  i  c r oe 1 e c t r ode  to 
reach  the  output  undistorted. 

The  partial  differential  equation  relating  voltage, 
position  and  time,  which  characterizes  the  general  distribut 
ed  R-C  line,  has  variable  coefficients  and  leads  to  a  situa¬ 
tion  which  is  easily  solved  only  in  a  small  number  of  cases. 
The  specialized  equation  for  the  particular  line  of  interest 
is  solved  and  the  steady-state  solutions  for  voltage  and 
current  presented.  On  this  basis  the  z-parameters  character 
izing  the  network  are  derived.  A  parallel  analysis,  which 
by-passes  the  solution  of  the  equation,  is  used  to  derive 
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the  voltage  transfer  function  of  the  network  which  yi 
information  regarding  the  time  domain  behavior  of  the 
and  forms  the  basis  for  a  stability  analysis. 
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CHAPTER  I 


I NTRODUCT I  ON 


1 .  1  Gene  r  a  1 

In  the  field  of  current  physiological  research,  and 
within  the  scope  of  biomedical  electronics,  a  topic  of  con¬ 
siderable  investigation  in  the  last  fifteen  years  has  been 

the  measurement  of  the  bioelectric  potentials  of  nerve  and 

(  1  -  9 ) 

muscle  cells.  The  measurement  of  these  potentials  has 

been  made  possible  by  the  development  of  the  glass  capillary 

(3) 

microelectrode  with  a  tip  diameter  of  0.5  microns  or  less. 

An  electrode  of  these  dimensions  allows  most  of  the  cells  of 
interest  to  be  penetrated  without  undue  damage  and  results 
in  highly  consistent  measurements. 

Measurements  with  small  electrodes,  however,  mean  com¬ 
plex  electronic  problems.  Subsequent  amplifying  stages  must 
meet  very  difficult  and  exacting  specifications^  ’  ’  ’  ^ 

and  the  distributed  R-C,  low-pass  character  of  the  micro¬ 
electrode  in  its  measuring  environment  introduces  distortion 
which  must  be  compensated.  With  the  application  of  new  cir¬ 
cuit  techniques  and  field-effect  transistors  with  low  noise 
figures,  present-day  physiological  amplifiers  are  as  suf¬ 
ficient  as  today's  state-of-the-art  pe rm i t s . ^ 1 1 ’ 1 ^ 

The  aspect  of  the  physiological  amplifier  which  con¬ 
cerns  this  particular  investigation  is  the  correction  of 
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the  waveform  distortion  which  is  introduced  by  the  micro¬ 
pipette  electrode.  Several  methods  are  available  to  the  de¬ 
signer  and  the  three  most  important  are  discussed  in  a  paper 
presented  by  Amatniek.^^  The  method  which  has  been  accept¬ 
ed  by  most  designers  and  which  has  met  the  specifications 
with  the  highest  degree  of  success  employs  input  capacity 
neutralization.  How ever,  certain  limitations  are  still  pre¬ 
sent  and  some  distortion  is  inherent  with  this  method.  With 
the  growing  importance  of  measuring  the  cellular  potentials 
with  increased  accuracy,  it  was  considered  within  reason  to 
search  for  some  improved  means  of  acheiving  the  desired  com¬ 
pensation. 

The  remainder  of  this  chapter  is  concerned  with  a 
more  detailed  analysis  of  the  problem  and  lays  the  basis  for 
a  scheme  of  compensation  which  employs  the  principles  of 
network  equalization. 

1 . 2  Physiological  Importance  of  Measurements 

All  living  tissue,  in  particular  the  tissues  of  nerve 
and  muscle  cells,  are  minute  generators  of  electricity.  Due 
to  the  differing  ionic  content  of  the  intracellular  and 
extracellular  fluids  a  potential  difference  (resting  poten¬ 
tial)  exists  across  the  walls  of  all  cells.  In  the  case  of 
excitable  tissue,  such  as  nerve  and  muscle,  action  potentials 
can  be  stimulated  and  are  associated  with  the  depolarization 
of  the  cell  wall.  These  transient  potentials  may  have  a  rate 
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of  rise  up  to  1000  v/sec^10,1^  (rise  time  of  100  ysec  for  a 

100  mv  pulse)  and  may  last  as  long  as  20  msec^)  with  1  msec 

being  typical.  The  action  potential  is  associated  with  the 
passage  of  an  impulse  along  the  nerve  fiber  or  with  the  con¬ 
traction  of  a  muscle  fiber,  and  forms  the  basis  of  much  re¬ 

search  in  neurophysiology.  The  measurement  of  electrical 
potential  gradients  across  living  membranes  has  aided  in  a 
better  understanding  of  the  ion  transport  mechanisms  involv¬ 
ed.  The  accurate  measurement  of  resting  and  action  poten¬ 
tials  constitutes  a  valuable  aid  towards  a  fuller  understand¬ 
ing  of  the  bioelectric  nature  of  living  tissue. 

1 . 3  The  Capillary  Microelectrode 

Both  resting  and  action  potentials  in  intracellular 
measurements  are  of  interest  in  cellular  research  today. 

The  capillary  microelectrode  is  used  widely  because  of  its 
ability  to  record  both  simultaneously.  However,  if  one  is 
interested  in  extracellular  records  and  rapid  changes  of  po¬ 
tential,  the  metal-tip  microelectrode  can  be  used  to  advant¬ 


age 
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Before  the  time  of  Ling  and  Gerard,  capillary  elec 

trodes  of  tip  diameter  2-5  microns  were  used  in  the  measure¬ 
ment  of  cell  membrane  potentials;  the  insertion  into  a  fiber 
caused  obvious  damage.  In  1949,  Ling  and  Gerard  developed 
and  perfected  methods  for  drawing  and  filling  microelectrodes 
of  well  under  one  micron  tip  diameter  and  were  rewarded  by 
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(b) 

obtaining  highly  consistent  membrane  potentials.  in  1953 

( 6 ) 

Alexander  and  Nastuk'  '  developed  a  mechanical  electrode 
puller  which  produced  a  more  consistent  group  of  electrodes 
at  improved  yield.  Microelectrodes  are  seldom  drawn  by  hand 
today . 

Capillary  microelectrodes  today  are  characterized  by 
the  following  properties: 

1.  Microelectrodes  are  usually  formed  from  a  hard 
borosilicate  glass,  usually  Pyrex  or  Phoenix,  to  en¬ 
sure  strength  and  minimum  chemical  reaction  with  the 
the  cellular  fluid. 

2.  A  two-stage  mechanical  (solenoid)  draw  is  used  to 
achieve  constant  and  repeatable  results. 

3.  The  ratio  of  lumen  to  wall  thickness  remains  app¬ 
roximately  constant  to  the  very  tip. 

b.  Capillary  has  a  fine  taper  (terminal  2mm)  with  a 
gradient  of  1/60  to  1/10. 

5.  Tip  outer  diameters  are  0.5  to  1  micron  to  ensure 
minimum  cell  damage. 

6.  Microelectrodes  are  filled  with  3  M  KC1  resulting 
in  comparatively  1 ow  tip  resistance  and  minimum  re 
action  with  the  cellular  fluid. 

7.  Electrode  resistance  may  vary  from  5  to  50  megohms, 
but  is  usually  around  20  megohms.. 

8.  Capacitance  of  the  electrode  is  approximately  1 
pf .  per  mm .  of  length. 


I 


* 

. 


' 

. 

a  a^mm  ^  kill  III 


(5) 


A  capillary  microelectrode  produced  from  a  two-stage 
mechanical  draw  is  represented  d i ag r amma t i ca 1 1 y  in  Fig.  1-1 


Stem 


Shoulder 


Shank 


Tip 


Pictorial  Representation  of  Microelectrode 

1 . k  Measuring  with  Microelectrodes 

To  study  the  electrical  activity  of  individual  cells 
the  fluid-filled  capillary  m i c roe  1 e c t rode  is  introduced  into 
the  interior  and  the  potential  difference  across  the  cell 
wall  determined.  The  tissue  is  bathed  in  Ringer  s  fluid 
with  the  bath  earthed.  A  metal  electrode,  usually  Ag-AgCl, 
is  placed  in  the  bath,  away  from  the  tissue,  where  the  chem¬ 
ical  environment  is  not  subject  to  variation.  To  compensate 
potentials  at  the  metal-fluid  junctions  another  Ag-AgCl  elec 
trode  is  connected  to  the  microelectrode  filled  with  3  M  KC 1 
The  use  of  3  M  KC 1  as  the  caoillary  electrolyte  ensures  low 
tip  resistance  and  a  low  diffusion  potential  at  the  fluid 
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tip  junction.  The  glass  capillary  mi croelectrode  serves  as 
an  e 1 e c t r o 1 y te  ' b r i d ge  between  the  interior  of  the  cell  and 
the  Ag-AgCl  electrode;  a  fluid  junction  is  established  at 
the  tip  within  the  cell.  Both  metal  electrodes  are  connect¬ 
ed  to  the  input  of  the  amplifier.  An  illustration  of  the 
measuring  setup  is  provided  in  Fig.  1-2.  The  electrode  is 
usually  mounted  as  close  as  possible  to  the  amplifier  input 
in  order  to  improve  the  high  frequency  response  of  the  re¬ 
cording  system.  Various  screening  schemes  are  employed  to 
reduce  stray  capacitances.  Excellent  references  concerning 

m i c roe  1 e c t rode s  and  their  use  are  Cater  and  Silver!^ 

(8)  (o) 

Gesteland,  et.  al,  and  Kennard. 


Figure  1-2 

Arrangement  for  Measuring  with  Microelectrodes 
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1  .  5  Previously  Proposed  Equivalent  Circuits 

To  facilitate  some  scheme  of  compensation  for  the 
distortion  that  the  low-pass  character  of  the  capillary 
m i c roe  1 ect rode  introduces,  the  designer  must  have  available 
some  equivalent  electrical  representation  of  the  electrode 
in  its  measuring  environment.  Two  researchers,  Woodbury^) 
and  Amatniek^^  have  presented  models. 

Woodbury's  model  assumes  that  the  resistance  of  the 
electrode  and  the  capacity  between  the  internal  and  external 
conducting  media  make  the  capillary  microelectrode,  electri¬ 
cally,  a  continuous  R - C  transmission  line.  The  model  is  pre¬ 
sented  in  Fig.  I  -  3  ( a )  and  the  numerical  values  are  based  on 
assuming  a  total  electrode  resistance  of  50  megohms  and  that 
the  electrode  is  conical  in  shape  with  a  slope  of  1/66. 
Capacitance  per  unit  length  is  constant  based  on  the  assump¬ 
tion  that  the  ratio  of  lumen  to  wall  thickness  remains  const¬ 
ant  throughout.  He  further  assumes  that  the  resistances  and 
capacitances  of  the  circuit  in  Fig.  I  —  3  ( a )  can  be  combined  to 
yield  the  more  apDroximate  circuit  of  Fig.  I  -3(b)  .  With 
stray  capacitances  in  the  range  of  2  ->  20  pf,  the  micro¬ 
electrode  is  characterized  with  a  time  constant  in  the  range 
of  0.2  ->  1  msec  and  provides  severe  distortion  to  an  action 
potential  with  a  rise  time  of  100  ysec. 
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Amatniek  presents  a  more  sophisticated  model  which  is 
derived  by  dividing  the  analysis  of  the  microelectrode  into 
three  segments:  region  (c),  the  portion  of  the  microelectrode 
within  the  cell,  contains  most  of  the  electrode  resistance; 
region  (d),  the  region  surrounded  by  external  fluid,  contains 
some  of  the  electrode  resistance;  and  region  (e),  the  portion 
in  air,  contains  very  little  electrode  resistance.  Distri¬ 
buted  capacitances  are  present  along  the  resistances  in  a 
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low-pass  manner.  The  capacity  and  resistance  of  region  (e) 
may  be  neglected.  Fig.  1-4  illustrates  Amatniek's  proposal. 


Tip-  Rc  R  d 


Figure  1-4 

Amatniek's  Equivalent  Circuit 

Previous  compensation  schemes  have  been  based  on  a 
simple,  two-element  R-C  low-pass  network  representation  sim¬ 
ilar  to  Woodbury's  as  illustrated  in  Fig.  I  -  3 ( b )  .  The  resis¬ 
tance  was  the  total  electrode  resistance  while  the  capacit¬ 
ance  was  a  lumped  sum  of  total  electrode  capacitance,  stray 
capacitance  and  amplifier  input  capacitance. 

1.6  M  e  t  h  ods  of  Reducing  Distortion 

The  action  potential  from  a  nerve  or  muscle  cell  may 
last  anywhere  from  500  ysec  to  20  msec  with  a  rise  time  as 
small  as  30  ysec.  Some  method  of  reducing  the  distortion 
introduced  by  the  low-pass  microelectrode  is  necessary  if 
the  high  frequency  changes  of  potential  are  to  be  followed 
accurately.  There  are  three  theoretically  feasible  methods 
for  restoring  the  undistorted  waveform:  decreasing  electrode 
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resistance,  decreasing  electrode  capacitance  (including  cap¬ 
acity  neutralization),  and  equalizing  the  overall  responseP^ 
1.6.1  Reduction  of  Capacitance  and  Resistance 

If  both  the  resistance  and  capacitance  of  the 
input  circuit  to  the  amplifier  can  be  reduced  suffi¬ 
ciently,  the  input  circuit  will  be  characterized  by  a 
fast  time  constant  (10  ysec  required  to  ensure  that 

the  recorded  action  potential  is  but  a  little  distort- 

(12) 

ed)  and  the  undistorted  waveform  can  be  restored. 

We  are  forced  to  tolerate  high  resistance  microelect¬ 
rodes  (in  the  order  of  20  megohms)  since  small  tip  di¬ 
ameters  and  low  shank  gradients  are  necessary  for  cell 
penetration  with  minimal  damage.  Capacitance  reduction 
has  been  accomplished  using  sophisticated  electrode 
fabrication  techniques,  shielding  of  input  leads  and 
special  amplifier  design  techniques.  By  far  the  most 
successful  means  of  reducing  input  capacity  has  been 
capacity  neutralization, 

1.6.2  Capacity  Neutralization 

Capacity  neutralization  is  based  on  the  follow¬ 
ing  simple  principle:  two  capacities  of  the  same  mag 
nitude  but  opposite  polarity,  when  paralleled,  will 
neutralize  or  cancel  each  other.  The  essential  ele 
ment  of  the  neutralized  amplifier  is  the  generation  of 
a  negative  capacity  at  its  input  terminals.  This  is 
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easily  obtained  via  the  Miller  effect  when  output-to- 
input  capacitive  feedback  is  used.  The  input  capacity 
is  easily  balanced  to  a  value  close  to  zero  by  adjust¬ 
ing  the  feedback  capacitor.  Continuous  compensation, 
which  may  be  necessitated  by  in-experiment  changes  in 
the  equivalent  input  circuit  (due  to  such  factors  as 
tip  breakage  or  tip  clogging),  is  easily  monitored. 
Restoration  with  capacity  neutralization  applied  has 

yielded  overall  time  constants  in  the  order  of  8  ysec 

(  8 } 

for  a  9  megohm  electrode.  ' 

However,  certain  limitations  are  still  present 
and  some  distortion  is  inherent  with  this  method. 

One  possible  cause  of  this  distortion  isthat  the  gen¬ 
erated  negative  capacity,  which  is  assumed  to  be  pure 
reactance  is  in  fact  a  complex  impedance  and  closely 
approximates  a  negative  capacity  over  only  a  limited 
bandwidth.  Also,  the  two-element  R-C  low-pass  equi¬ 
valent  circuit  is  a  very  simplified  and  approximate 
version  of  the  actual  network.  Even  if  the  input  cap¬ 
acity  of  the  amplifier  were  in  fact  neutralized,  there 
still  could  be  waveform  distortion  due  to  the  remain¬ 
ing  distributed  constants.  In  this  case  the  best  neg¬ 
ative  capacity  adjustment  might  be  different  from  an 
exact  neutralization  of  input  capacity. 
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.6.3  Equal  ization 

Another  method  of  compensation,  discussed  by 
Amatniek^1)  and  used  by  Woodbury  is  that  of  net¬ 

work  equalization.  This  method  involves  the  introduc¬ 
tion  into  the  amplifier  of  a  transfer  characteristic 
inverse  to  that  of  the  low-pass  capillary  microelect¬ 
rode  producing  a  flat  overall  response.  In  Fig. 1-5 
a  block  diagram  of  an  equalization  circuit  is  shown 
where  the  reciprocal  transfer  ch a r a c te r  i  s t  i  c  of  the 
m i c r oe 1 e c t ro de  is  introduced  between  two  stages  of  the 
amplifier.  Previous  usage  of  this  method  assumed  the 
two-element  equivalent  circuit  as  shown  and  V/oodbury 
was  able  to  realize  overall  system  time  constants  in 
the  order  of  Ad  ysec  for  a  22  megohm  electrode. 


Figure  1-5 

Compensation  by  Equalization 

The  generation  of  the  reciprocal  transfer  function 
can  be  accomplished  by  placing  a  circuit  with  a 
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transfer  characteristic  approximately  equal  to  that 
of  the  m i c roe  1 e c t rode  in  an  i  n p u t  -  to - ou t p u t  feedback 
loop  of  an  operational  amplifier.  Even  with  the  op¬ 
erational  realization  of  such  a  circuit  the  method 
suffers  one  major  drawback.  If  the  electrode  resis¬ 
tance  changes  due  to  clogging,  tiD  breakage,  or  chang¬ 
ing  electrolyte  concentration  the  low-pass  cutoff  fre¬ 
quency  also  changes,  and  therefore  the  equalizing  net¬ 
work  must  be  readjusted.  If  the  equalization  is  not 
continuously  monitored  and  the  change  takes  place  in 
the  middle  of  an  experiment,  incorrect  results  may  be 
recorded  unknowingly. 

1 . 7  Proposed  Scheme  of  Waveform  Restoration 

Since  the  fluid-filled  capillary  microelectrode  be¬ 
haves  electrically  like  a  continuous  R-C  transmission  line 
(section  1.5),  it  is  considered  feasible  to  achieve  wave¬ 
form  restoration  by  employing  the  principles  of  network 
equalization  using  thin-film  distributed  parameter  circuits. 
If  feasible,  such  an  evolution  would  have  the  advantage  of 
producing,  with  more  accuracy,  the  inverse  transfer  charac¬ 
teristics  required  to  equalize  the  overall  response.  Also, 
equalization  is  superior  to  capacity  neutralization  from  the 
viewpoint  of  noise  considerations.^^  The  proposed  equaliza¬ 
tion  scheme  is  illustrated  in  Fig.  1-6.  Fig.  I  -6(a)  illus¬ 
trates  the  entire  scheme  in  block  form.  corresponds  to 
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Figure  1-6 

Proposed  Scheme  of  Waveform  Restoration 

the  generated  cell  potential  and  V is  the  restored,  undis¬ 
torted  waveform  amplified  by  the  factor  A-|A2.  V ^  can  be  fed 
directly  to  an  oscilloscope  or  may  undergo  further  amplifica¬ 
tion  before  display.  Fig.  I  -  6  ( b )  illustrates  h  ow  the  block 
1/H(s)  is  acheived.  The  required  condition  is  contained  in 

the  following  equation: 

A  1  H  (  s ) V  1  -  H(s)V4  =  0 

This  condition  can  be  met  with  the  use  of  an  operational 
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amplifier  characterized  by  high  input  impedance  and  high  dif¬ 
ferential  voltage  gain.  The  output  of  the  operational  ampli¬ 
fier  will  be  an  amplified  version  of  V  ^  and  the  undistorted 
waveform  will  be  restored. 

1 . 8  Subsequent  Analysis 

In  the  analysis  to  follow,  a  proposal  of  an  equaliza¬ 
tion  scheme  for  compensation  is  presented  and  its  feasibility 
tested.  As  part  of  the  initial  development,  existing  models 
of  the  microelectrode  are  presented  and  used  as  a  basis  for 
establishing  a  convenient  and  accurate  electrical  equivalent 
of  the  probe.  The  majority  of  this  paper  is  concerned  with 
the  mathematical  description  of  the  electrode,  both  from  the 
point  of  view  of  its  two-port  parameters  and  its  voltage 
transfer  function  H(s).  The  latter  Is  useful  in  establishing 
the  stability  of  the  proposed  equalization  scheme. 
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CHAPTER  I  I 


THE  EQU I  VALENT  CIRCUIT 


2 . 1  General 

The  mathematical  analysis  of  any  physical,  electrical 
system  requires  the  choice  of  some  equivalent  electrical 
representation  of  the  system.  The  equivalent  circuit  is 


usually  chosen  to  be  as  simple  as 

the  1  i 

m  i  t  s 

o  f 

accu  racy 

re  - 

q  u i  r e  d  will 

allow.  This  ensures, 

h  o  pe  f  u 

Hy, 

a 

possible, 

i  f 

not  simple, 

mathematical  analysis 

of  the 

s  y s  tern 

yielding 

the 

desired  design  information. 

The  system  of  interest  in  this  investigation  is  the 
e 1 e c t ro 1 y te - f i 1 1 e d  capillary  microelectrode  in  its  measuring 
environment.  In  all  previous  work,  designers  have  relied 
on  the  two-element  R-C  low-pass  network  (Figure  I  -  3  b )  as 
this  system's  equivalent  circuit.  All  forms  of  reducing 
distortion  have  been  based  on  this  simplified  representation. 
Closer  to  the  problem  at  hand,  Woodbury  used  network  equaliza¬ 
tion  based  on  this  simple  model  to  approximately  restore  the 
undistorted  waveform. ^  ^  He  was  quite  successful  with  slow 

rises  of  potentials  but  fast  rising  potentials  reached  the 
display  distorted. 

The  selection  of  a  model  in  this  investigation  is 

based  on  the  foil o w ing  criteria: 

1.  A  more  accurate  representation  of  the  system  than 

was  assumed  by  previous  investigators  with  regards 
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to  waveform  restoration. 

2.  Possibility  of  fabrication  using  thin-film  dist¬ 
ributed  parameter  networks. 

3.  Fabricated  circuit  of  the  model  should  lend  itself 
to  a  method  of  continuous  equalization  should  the 
characteristics  of  the  physical  system  change. 

The  remainder  of  this  chapter  is  devoted  to  the 
adoption  of  an  improved  model  and  a  mathematical  representa¬ 
tion  of  its  distributed  parameters. 

2 . 2  Adopted  Equivalent  Circuit 

From  the  considerations  of  Section  2-1  the  following 
equivalent  circuit  has  been  adopted  as  a  more  accurate 
representation  of  the  microelectrode  in  its  measuring  en- 
v i r on  me  n  t : 


r(x)  is  resistance  variation  along  the  electrode 
c(x)  is  capacitance  variation  along  the  electrode 
C.  is  the  amplifier  input  capacitance 


Figure  11-1 

Adopted  Equivalent  Circuit 
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This  model  is  based  on  Woodbury's  consideration  that  the 
electrode  is  equivalent  electrically  to  a  continuous  R-C 
transmission  line.  C.  is  included  to  account  for  the  cap¬ 
acitance  contributed  to  the  system  by  the  amplifier  and  the 
interconnecting  leads. 

This  equivalent  circuit  is  definitely  superior  to 
the  two-element  model,  but  essentially  consistent  with 
Amatniek's  proposal.  Due  to  its  simplicity  and  contin¬ 


uity,  the  present  model  should  lend  itself  to  thin-film 
fabrication  and  some  form  of  continuously  monitored 
equal ization. 

2-3  Electrode  Capacitance  Variation 

The  determination  of  the  capacitance  variation 
along  the  electrode  is  based  on  the  following  known  facts: 

1.  The  depth  of  Ringer's  fluid  above  the  tissue  is 
about  2  mm. 

2.  Glass  is  usually  either  Pyrex  or  Phoenix. 

3.  Lumen  is  normally  2/3  of  outer  diameter. 

4.  Ratio  of  lumen  to  wall  thickness  remains  constant 
throughout  with  the  ratio  varying  less  than  5%. 

At  any  particular  point  of  the  microelectrode  immersed  in 
Ringer's  fluid,  the  electrode  has  the  capacitive  properties 
of  a  coaxial  cylinder.  This  is  illustrated  in  Fig.  11-2. 
The  coaxial  cylinder  is  characterized  by  the  capacitive 
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Figure  11-2 

Electrode  Cross-Section 


distribution  c ( x) =2 1 n ( r/a  )  ,  where  e  is  the  dielectric 
constant  of  the  insulating  medium.  With  the  slope  of  the 
capillary  gradual  and  the  ratio  of  lumen  to  outer  diameter 
constant  throughout  the  length,  the  equation  for  the  cap¬ 
acitance  of  a  coaxial  cylinder  is  applicable  to  the  electrode. 
This  indicates  that  the  distributed  capacitance  of  the  elec¬ 
trode  will  be  constant  (per  unit  length).  This  is  in  ac¬ 
cordance  with  the  assumptions  of  previous  investigators. 

2  .  4  Electrode  Resistance  Variation 

In  the  determination  of  the  resistance  variation 
along  the  electrode  the  following  assumptions  are  made: 


1.  All  electrode  resistance  is  due  to  the  thin  column 

of  electrolyte  within  the  capillary. 

2.  The  electrode  is  conical  in  form  with  constant 

slope  in  the  shank  region. 

3.  Electrolyte  resistivity  is  constant  throughout. 


These  assumptions  provide  the  basis  for  a  mathematical 
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description  of  the  electrode  resistance.  The  accuracy  of 
this  treatment  will  be  sufficient  within  the  limits  of  the 
assumed  model. 

The  derivation  of  the  per  unit  length  resistance 

r(x)=iL!i  of  the  conical  column  is  based  on  the  well  known 

dX  L 

relationship  governing  the  resistance  of  a  conductor:  R=p^\ 

A  geometrical  representation  of  the  column  is  given 

in  Figure  I  I  -  3  a  •  Aside  from  the  obvious  area  relationships 


A 


Figure  11-3 

Geometrical  Representation  of  Electrolyte  Column 
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s  h  ow  n  in  the  diagram,  the  following  equations  may  be 

radius  of  column  at  x 


a (x)  =  aQ  +  ax 


A(x)  =  ir(a0  +  ax)  cross-sectional  area  at 
A  differential  element  of  the  column  is  presented  in 
Figure  I  I  -  3  b  .  For  the  element: 


r  (  x )  =  iL!i  =  1  i  m  AR_ 
d*  Ax->0  Ax 


written: 


x 


It  can  easily  be  shown  that: 

r  (  x )  =  -  =  - - 

A  (  x)  tt  (  aQ+ax)  z 

(for  details  consult  Appendix  I) 

The  function  r(x)  above  describes  the  resistance  variation 
of  the  electrolytic  column  within  the  capillary  of  the 
m i c r oe 1 e c t r o de .  The  space  parameter  x  is  always  positive 
and  increases  along  the  axial  length  of  the  microelectrode, 
away  from  the  tip. 

2 . 5  S  umma  ry 

In  this  chapter  an  equivalent  circuit  of  the  elect¬ 
rolyte-filled  capillary  microelectrode  in  its  measuring 
environment  has  been  adopted,  based  largely  on  the  work  of 
previous  investigators.  The  model  consists  of  a  R-C  low-pass 
distributed  parameter  network,  representing  the  microelec¬ 
trode,  shunted  by  a  single  capacitor  representing  all  cap¬ 
acitive  contributions  exterior  to  the  electrode.  The  resis¬ 
tance  and  capacitance  variations  of  the  distributed  parameter 
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network  have  been  derived  and  are  presented  as  follows: 


c  ( x) 
r  (x) 


c 


o 


P 


it  (  a  +ax)  2 
o 


{2-1} 


{2-2} 


CHAPTER  I  I  I 


MATHEMAT I  CAL  ANALYS I S 


3  .  1  General 

In  the  analysis  to  follow,  the  equivalent  circuit  of 
the  microelectrode  in  its  measuring  environment  is  divided 

into  two  functional  blocks.  This  is  illustrated  in  Figure 
lll-l.  This  division  is  convenient  for  two  reasons: 

1.  It  isolates  and  simplifies  the  analysis. 

2.  It  constitutes  a  logical  division  with  respect  to 
some  ensuing  scheme  of  continuous  equalization,  since 
the  elements  of  each  block  are  subject  to  independent 
variation. 

Due  to  the  complexity  inherent  to  the  analysis  of  distri¬ 
buted  parameter  networks,  the  majority  of  this  chapter  is 
concerned  with  the  analysis  of  the  first  functional  block. 
The  consideration  of  the  two  blocks  cascaded  results  in  a 
description  of  the  complete  network. 
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Figure  lll-l 

Microelectrode  Equivalent  Circuit 
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In  this  investigation  the  goals  of  the  analysis  are 
two- f o 1 d : 

1.  To  provide  the  basis  upon  which  the  feasibility 
of  the  compensation  proposal  can  be  tested. 

2.  To  provide  information  which  may  be  of  interest 
to  investigators  concerned  with  the  properties  of 
non-uniform  distributed  R-C  lines. 

In  the  analysis  to  follow  the  characterizing  partial  dif¬ 
ferential  equation  for  voltage  along  the  line  is  solved. 

The  resulting  product  solution  is  used  to  resolve  the  sinu¬ 
soidal  steady-state  solutions  for  both  line  voltage  and 
current.  These  solutions  are  employed  to  derive  the  z-para- 
meters  of  the  line.  The  derived  z-parameter  matrix  contains 
all  the  information  necessary  towards  the  derivation  of  any 
two-port  network  parameter  of  interest.  A  recently  evolved 
method  of  analysis^30)  is  used  to  derive  the  network  func¬ 
tion  of  the  system,  where  H  ( s )  is  expressed  in  terms  of  a 
uniformly  convergent  series  of  s.  The  coefficients  of  this 
series  are  used  to  extract  information  regarding  the  line's 
time  domain  behavior.  The  network  function  H(s)  will  form 
the  basis  for  subsequent  tests  which  will  be  necessary  in 
determining  the  stability  of  the  proposed  compensation 


scheme. 
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'3*2  Analysis  of  the  Distributed  R-C  Line 


Part  I  -  The  Differential  Equation 


The  first  functional  block  of  the  complete  network 


represents  a  tapered  distributed  R-C  line  with  the  following 
resistance  and  capacitance  distributions: 


c(x)  =  cQ 


{3-1} 


(3-2} 


With  the  increased  interest  in  microelectronics  and  integ¬ 
rated  circuits,  the  analysis  of  distributed  parameter  net- 

(14-31 

works  has  received  the  attention  of  many  investigators 
in  recent  years. 

In  the  analysis  to  follow,  the  traditional  approach 
will  be  employed.  The  analysis  will  be  based  upon  the  solu¬ 
tion  of  the  line's  ch a r a c te r i z  i  n g  differential  equation. 

From  the  complete  solution  the  sinusoidal  steady-state  solu¬ 
tion  is  readily  derived.  The  steady-state  solution  is  in 
turn  used  to  derive  any  required  set  of  the  network's  two- 
port  parameters.  In  this  case,  the  network's  z-parameters 

are  presented. 

3-2.1  The  Differential  Equation 


Kaufman  and  Garrett^2^  have  shown  that,  if 
attention  is  restricted  to  one  dimensional  current 
f 1 ow ,  the  general  nonuniform  distributed  R-C  line  is 
characterized  by  a  second-order  partial  differential 
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equation  with  variable  coefficients.  The  derivation 
is  presented  in  Appendix  II  with  the  following  equa¬ 
tions: 


9  2  v 
9  x2 


— L.  lidisi  i*  =  r(x) 

r ( x)  dx  9x 


{3-3) 


9  2  i _ 1_  d  c  (  x )  9_i_ 

9x2  c(x)  dx  9x 


r  (  x)  c  (  x) 


Ai 

9  t 


(3-4} 


Both  voltage  (v)  and  current  (i)  are  functions  of 
both  time  (t)  and  distance  along  the  line  (x).  The 
solution  of  these  equations  in  closed  form  for  arbi¬ 
trary  r(x)  and  c(x)  has  not  been  found;  even  the 

(25) 

steady-state  solution  is  too  complicated.  A 

straight  forward  solution  exists  for  the  uniform  line 
(no  taper)  while  other  solutions,  if  possible,  are 
quite  involved.  The  special  case,  where  the  product 
r(x)c(x)  is  a  constant,  has  received  considerable 
treatment  because  of  the  ease  of  fabricating  such  a 
line. 

In  this  investigation  the  solution  of  { 3  ~  3 } 
with  the  resistance  and  capacitance  distributions  of 
(3-1}  and  (3-2}  is  of  interest.  With  the  substitu¬ 
tion  of  (3"l}  and  (3"2}  : 

_ 1_  d  r (x)  _  2a  { 3 -  5 } 

r  ( x)  dx  aQ+ax 
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r(x)c(x)  =  - — 

7T  ( aQ+ax)  2 

Therefore  { 3 - 3  1  becomes: 

9 2  v  2  3v_  =  _ pco  3_v 

3x2  aQ+ax  3x  Tr(a0+ax)2  3t 

This  can  be  written  in  the  form: 


(3-6} 


(3-7) 


IT 


P  C 


(a0+ax)  - 


3  2  v  (  x  ,  t )  2  7Ta 


Pc, 


( aQ+ax) 


3  v  (  x  ,  t )  _  3  v  ( x  ,  t ) 

3x  3 1 


{3-8} 


3.2.2  Solution 

A  solution  will  be  attempted  of  the  form: 

v ( x  ,  t )  =  T ( t )  X  (  x )  =  V ( t )  V ( x)  {3-9} 

Substitution  of  the  trial  solution  into  the  equation 
( {3- 8})  and  separation  of  variables  results  in  two 
separate  ordinary  differential  equations,  one  with 
time  t  as  the  independent  variable  and  the  other  with 
distance  x  as  the  independent  variable: 


TT 

pco 


( a  Q+ax) 2 


1  ill 

X  dx2 


2  7Ta 

Pco 


(a0+ax) 


i 

X  dx 


k2 


(3-10} 


1  H  =  -  k2  {3-11} 

T  dt 


This  separation  into  two  ordinary  differential  equations 
indicates  that  a  product  solution  does  in  fact  exist. 

The  parameters  k2  which  satisfy  the  equations  and  as¬ 
sociated  boundary  conditions  are  termed  eigenvalues 
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and  the  resulting  functions  X  ( x )  and  T(t)  are  termed 
eigenfunctions.  The  final  solution  to  an  eigenvalue 
problem  may  be  expressed  in  series  form  composed  of 
an  infinite  number  of  normal  modes . ^3) 
the  solution  T ( t ) 

The  equation  may  be  expressed  as: 

”  +  k 2  T  =  0  {3-12} 


The  solution  is  easily  resolved  as: 

-k2 1  + 

V ( t)  =  T(t)  =  A  e  ;  A  =  T(t=0  )  (3-13) 

the  solution  X(x) 

The  equation  may  be  expressed  as: 
d2X  dX  pcQ 

(a  +ax)  2  - -  +  2a(a  +ax)  —  +  -  k2X  =  0  (3-1*0 

dx2  dx  7T 


This  equation  has  the  form  of  Legendre's  Linear 
Equation  and  can  be  reduced  to  an  ordinary  differen¬ 
tial  equation  with  constant  coefficients  by  using  the 

r  ,  1  .  (32) 

following  transformation  : 

ez  =  (a  +ax) 
o 

The  details  of  the  transformation  are  presented  in 
Appendix  III.  In  the  transformed  space  the  equation 
b e come s  : 


d 2  X  dX  pcQ 

-  +  -  +  -  k 2  X  =  0  ;  X  (  z )  «~>X  (  x)  (3-15) 

dz2  dz  Trot 2 
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Y  z 

Assuming  a  solution  of  the  form  X=e  yields  the 
following  characteristic  equation: 

y2  +  y  +  —  k2  =0  {3~16} 

TTOt 2 

Solution  of  this  equation  yields  the  characteristic 
roots : 

y  =  i  ( - 1  ± G  1  )  {3-17> 

1  2 

pc  1 

0'  =  0  -  i.k2-4)2  {3-18} 

TTOt 


The  solution  in  the  z-space  takes  the  form: 


X(z) 


B  e 


Y  1  z 


+  C  e 


Y?z 


(3-19) 


where  B  and  C  are  constants  dependent  on  the  boundary 
condi  t i on  s . 

The  solution  can  also  be  expressed  as: 


X  (  z) 


-iz  ,  D  iG'z  A  r  -iG'z  v 
:  -  (  B  e  +  C  e  ) 


{3-20} 


Consider  the  identities 


y  G  1  z 


cos 


h  4 6 ' z  +  sinh  i G 1 z 


-4G  1  z 


cosh  iG 1 z  -  sinh  20'z 


Using  these  relationships,  equation  {  3 _  2  0 }  becomes 


.  • 


(30) 


where:  c  =  B  -  C  ;  c  =B+C 

1  2 

This  form  of  the  solution  will  prove  more  convenient 
in  the  subsequent  analysis.  Since  eZ  =  aQ  +  ctx,  the 
solution  in  the  x-space  is: 

V(x)  =  X(x)  =  (aQ  +  ax)  2  [  c^sinh(j)(x)  +  c^cosh^Cx)  ] 


where: 


cj>  (  x)  =  10  1  1  n  (  aQ+ax) 

0 1  =  (1  -  4k2^-)^ 

2 

Trot 

{3-21 } 


Substituting  (3“13}  and  (3~21}  into  (3-9)  yields  a 
solution  (3-22): 


e  -  k2  t 

v  (  x  ,  t )  =  - - -  [c^s  i  n  h  4>  (  x)  +  c^cosh^x)] 

( a  +ax) 2 


(3-22} 


where:  c„  =  A  c,  ;  c  -  A  c 

3  1  A2 

the  complete  solution 

A  solution  of  the  line's  characterizing  differential 
equation  ({3-8})  has  been  found  and  is  expressed  in 
{3-22}.  The  general  solution  may  be  expressed  as: 


v(x,t)  =  c  +  £  X(k.,x)  T(k.,t) 

(J  j  *  ' 

Therefore  the  complete  solution  to  {3“8}  may  be  ex¬ 
pressed  as: 

.  -k?t 

=  C  +  l — - - '■ - -  [c_s  i  nhtj)  (  x)  +  c^cos  h<j>  (  x)]  {  3  -  2  3  > 

®  i  ( a  +ax) 2 
o 


v(x, t) 


' 


. 


’ 


. 


(31) 


where  c ^ ,  c ^  and  c  ^  are  constants  determined  by  the 
conditions  to  which  the  line  is  subjected  and  the 
values  kj  represent  the  eigenvalues. 

3.2.3  Sinusoidal  Steady-State  Solution 
the  voltage  solution 

For  the  sinusoidal  steady-state  response: 


v  (  x  ,  t ) 


Re 


[c  +  V(x)eSt] 
l  0  J 


s  =  j  a) 


From  (3-21}  : 

V(x)  =  (a  +ax)  2  I  c.s  inh(j)(x)  +  c0  cos  h(|)  ( x) 

O  I  Z 

From  {  3 -  2  3  )  recognize:  s  =  -  k  2  =  j  co 

Therefore:  k2  =  -joo  {  3  ~  2  4  } 

1 

The  sinusoidal  steady -state  solution  of  line  voltage 
is  recognized  as  the  space  solution  { 3 —  2 1  }  held  to 
the  confines  of  {  3  "  2  4  }  .  Therefore,  in  the  sinusoidal 
steady-state: 


v 


(x  ,  t) 


_  1 


s  s  s 


=  V  ( x)  =  (aQ+ax)  2[Cls 


n  h<j>  (  x)  + 


c^  cos  h  4>  ( x )  ] 


where: 


(p  ( x) 

e 1 


id  1  1  n  (  a 
(  1  +  j  co 


+ax) 

o 

!tPco) 

TTOt  2 


JL 

2 


{3-25} 


' 


(32) 


the  current  solution : 

It  has  been  shown  (Appendix  II)  that  the  current  in 
the  line  is  related  to  voltage  by: 


dV  ( x ) 

d  x 


1  3  v 


r ( x)  3  x 


i  =  i  (  x ,  t) 
v  =  v  (  x  ,  t) 


For  the  sinusoidal  steady-state  response: 

i(x,t)  =  Re  [l(x)est]  ;  s  =  joo 
The  refore : 

1  dV ( x) 


I  (x) 


r ( x)  d  x 


Evaluating  —  V ( x) 
dx 


d  V  ( x) 


d  x 


[ 


ci 


dx  (aQ+ax) 


— s  i  n  h  ch  (  x )  + 


c2 


{3-26} 


(3-27) 


(a  +ax) * 


cos  h  <j)  (  x) 


(3-28} 

Expanding  the  differentiation  yields: 

-(a  +ax)  ”  2  F  (  c  0  '  -  c  , )  s  i  nh<J>  (x)  +  (  c  <  0  '  -  c  )  cos  h  <J>  (  x )] 

o  o  L  ?  1  '2 


i  (x  ,  t) 


Therefore,  in  the  sinusoidal  steady-state: 

=  -  ™(a  +ax)?[(c0'-c1)  s  i  n  h(f)  (  x)  +  (  c  .  0  '  -  c  )  cos  h<j>  (  x)] 

s  s  s 


where:  cf>(x)  =  i0  '  1  n  (  aQ+ax) 

0'  =  (1  +  joj— -f)2 

ira 

Equations  (3-25)  and  { 3 “ 2  9  >  represent  the  si 
steady-state  solutions  for  the  line.  These 
form  the  basis  for  the  following  development 
the  impedance  matrix  for  the  network. 


(3-29) 

n  us  o  i  d a  1 
equations 
deriving 


1 


m 

■ 


(33) 


3.2.4  Network  z  Parameters 

In  section  3.2.3  the  following  expressions  were 
derived  for  the  voltage  and  current  sinusoidal  steady- 
state  responses: 

V  ( x)  =  (  a  Q+ax)  2  [c  sinh<|>(x)  +  c2  cos  h  <p  (  x)]  {  3  -  3  0  } 

I  (  x )  =  -  — —  (  a  +ax)  2  [  (  c9  0  '  -  c  .  )  s  i  n  h  cj)  (  x).  +  (  c  0  1  -  c  )  cos  h<|>  (  x)] 

2  p  o  z  i  i  z 

(3-31) 

Consider  the  following  two-port  network  with  the 
boundary  conditions  as  shown: 


1 


1 

» 

1 

— 6 

1 

1 

1 

1 

\ 

O - 

1 

1 

1 

1 

t 

1 

1 

1 

1 

1 

NETWORK 

t 

1 

1 

1 

1 

1 

1 

1 

0 _ 

- O 

1 

1 

t 

1 

1 

x=0  x=A 


Figure  I  I  I  -  2 
Two  Port  Network 

With  the  voltage  and  current  equations  governing  the 
network  known,  the  z-parameter  matrix  can  be  formed: 


o  r 


— 

V 

z 

z 

1 

1  1 

1  2 

V 

z 

z 

2 

2  1 

2  2_ 

-  [Z][l 


{3-32} 


V] 


{3-33} 


- 


' 


> 


. 


(3*0 


Consider  the  boundary  conditions: 

at  x=0  V(x)=V1  at  x  =  X  V(x)=V2 

I  ( x )  =  I  1  -  I  (x)  =  I  2 

Substituting  these  conditions  into  { 3  ~  3  0  }  yields: 

V.|  =  V(0)  =  a  2[c^'sinhcj)(0)  +  c2cosh(j)(0)3  (3-34} 

V2  =  V  (  X)  =  (aQ+ax)  2  [c  ^  s  i  n  h(J)  (  X )  +  cos  h  cj>  (  X)  ]  (3_35) 

where: 

({)  (  0  )  =  i0  1  1  naQ 

<J)  (  X)  =  \ 0  1  1  n  (  aQ+ax) 


Expressing  { 3  ~  3  )  and  {3“35)  in  matrix  form: 


d 

d 

1  1 

1  2 

d 

d_ 

2  1 

22 

In  matrix  notation: 

v]  =  [d]  [c 


where: 

d  .  ,  =  a  2  sinh  ( 0 ) 

1  1  o 

d10  =  a  ~  2  cosh  (j)(0) 

1  £■  O 

d  =  (a  +aX)  2  s  i  n  h  4>  ( A ) 

2  1  o 

d  =  ( a  +aX)  2  cos  he})  ( X) 

2  2  o 


{3-36} 


(3-37) 


(3-38) 


■ 


(35) 


Substituting  the  boundary  conditions  into  { 3  -  3  1  } 
yields: 


I  =  1(0)  =  -  ^-a  2  [(  c  0  '  -  c  )  s  i  nh(j)  (  0  ) 

1  zp  o  ^  1 

+  (  c  1  0  '  -  c2  )  cos  h(j)  (  0  )]  {3“39} 

I  =  I  (A.)  =  ^L(a  +  aA)2[(co0'-c1)sinh<})(A) 

2  2  p  o  2  1 

+  (  c  ^  0  c  )  cosh(j>  (  A)]  {  3 ~  ^ 0  } 


Rearranging  {  3  “  3  9  3"  and  {  3  “  0  }  yields 


T TCZ  \ 

- — a 
2  p  o 


[sinh<j>(0)  -  0  1  cos  hcj)  (  0  )  ]  c  ^ 

+  [coshcj)(0)  -  0  !  s  i  n  he})  (  0  )]  c^ 
[sinh<j)(A)  -  0  '  cos  h  (})  (  A )]  c  ^ 

+  [coshc})(A)  -  0  1  s  i  n  h  <J)  (  A )] 


3^(a  +aA) 
2  p  o 


Expressing  { 3  ~  ^ 1 >  and  { 3 “ ^  2 }  in  matrix  form: 


(3-41} 


(3-^2} 


a 

a 

1  1 

1  2 

a 

a 

2  1 

22 

where 


a 

IT  Ct 

1 

2 

■  a 

[s 

i  n  h  cj)  (  0  )  - 

0  1  cos 

1  1 

2p 

0 

a  1  2 

7T  Ot 

2p 

i 

■  a5 
0 

1 — 1 
0 

0 

l/> 

z r 

-©- 

0 

0  '  s  i  n 

a  „  , 

=  _ 

ira 

( a 

+  aA  )  2  [s  i  n 

h  cj>  ( A ) 

2  1 

2  P 

0 

a  „ 

—  - 

7TCX 

( a 

+  aA )  2  [cos 

h  4>  ( A ) 

22 

2  P 

0 

(3-43) 


(3-H) 


In  matrix  notation: 


I]  =  [A][c 


o-'ts} 


♦ 


. 

. 


(36) 


Consider  the  following  matrix  development: 


n  ow 

V] 

V] 

1] 

[A]-’  1] 

11  11  11  11 

1 — 1  1 — 1  1  1  1  1 

L>|  ^3=^  |Q[  |M| 

1  1 - 1  1 - 1  1 - 1 

— *  O  O  — 

s 

1 — 1 

1  e  t 

1 - 1 

CD 

1 

1 — 1 
<l 

11 

therefore 

c] 

-  [B][l 

and 

V] 

■  [D][B][l 

the  re  f 0  re 

1 - 1 

|N| 

=  [D]  [8] 

(3-^6) 


The  next  step  in  the  development  will  be  the  formation 

of  matrix  [b]  which  is  the  inverse  of  matrix  [a]  . 

(33) 

For  the  2  x  2  array:7 


1 

de  t  [a] 


(3-47) 


where: 


d e  t  [a]  =  a  a  -  a  a 

1122  1221 


Evaluating  det 


a 


a 

1  1  22 


a  +  aA)l 
o 


+ 


1 

2 

0  ' 

0  1 

0  ' 


sinh4>(0)coshc{>(A) 

2sinh<H'Ocosh(j)(0) 

cosh  (f)  (  0  )  cos  h(j>  (  X) 
s  i  nh<j>  (0)  s  i  nh<j>(  A) 


(37) 


a 


a 

12  21 


[a  ( a 
u  o  o 


+  aX )  ] 


cosh(J)(0)sinh(})(X) 


+  0l2sinh4>(O)coshcj)(X) 


-  0'  cos  h  cj)  (  0  )  cos  h  <j)  (  X ) 

-  0'  sinh4>(0)sinh<})(X) 

det[A]  =  -  (|^)2  [a0(a0+aX)]i  (1  -  0'2) . 

•  •  •  •  [s  i  n  h  <j)  (  0 )  cos  h  cj)  (  X )  -  cos  hcf)  (0)  s  i  n  h  4>  (X)]  {  3“  ^8} 


Using  the  identity: 

sinh(x  -  y)  =  sinh  x  cosh  y  -  cosh  x  sinh  y 
Equation.  {  3  -  ^  8  }  becomes: 

de  t  [A]  =  -  (j~)  2  [aQ(ao+aX)]  2  (  1  -0  1  2)  s  i  nh  [(f)(0)  -  4>  (  X )  ] 

{ 3- 49  > 

where: 

ao  x 

<K0)  -  4>(X)  = 


Note:  From  Appendix  IV 

det[Aj  =  0  if,  and  only  ?  f  oo  =  0 
therefore  [b]  is  defined  for  all  conditions 

excep  t  a)  =  0  . 


. 


(38) 


Using  (3"^},  (3-^7)  and  (3"^9)  yields  the  parameters 

of  matrix  [b]  : 

_  a22  _  2p  coshcj)(X)  -  0'sinh(j>(X) 

^  d  e  t  [a]  ttcx 


1  2 


-a  i  £  -2p 

Tra 


det [A] 


a^(1-0'  2)sinh[4>(0)  -4>(X)] 
o 

0'sinhcJ>(O)  -  coshc{>(0) 

(a  +aX )-( 1  -  0  '  2 )  s  i  n  h  [$  (  0  )  -  4>  ( X )  ] 


-  a 


21  2p  0'cosh(j)(X)  -  sinhcj)(X) 


2  1 


det  [a] 


7T  Ot  X 


1  1  “2p 


a^(  1-0'  2)  s  inh  [cf>(0)  -  (^(X)] 

s  i  nhc(>  (0)  -  0  '  cos  hcj)  (0) 


22  det[A]  Tta  (a  +aX )  2  (  1  -  0  1  2  )  s  i  n  h  [(J)  (  0  )  -  <j)  ( X )  ] 


Now  recall: 


[z]  =  [d]  [b] 


{3-50} 

(3-51) 


o  r  : 


z 

z 

d 

d 

b 

b 

1 1 

1  2 

1  1 

1  2 

1  1 

1  2 

z 

z 

d 

d 

b 

b 

2  1 

22 

2  1 

22 

2  1 

22 

Expanding  the  product  yields 


=  d 

b 

+ 

d 

b 

1  1 

1  1 

1  1 

12  21 

=  d 

b 

+ 

d 

b 

1  2 

1  1 

1  2 

12  22 

=  d 

b 

+ 

d 

b 

2  1 

2  1 

1  1 

22  2  1 

=  d 

b 

+ 

d 

„  b 

22 

2  1 

1  2 

22  22 

(3-52} 


(3-53) 


Using  the  identities: 


sinh(x  -  y)  =  s i nh  x  cosh  y  -  cosh  x  sinh  y 

sinh  x  sinh  y 


cosh (x  ~  y)  -  cosh  x  cosh  y 


•» 


- 

and  some  algebraic  manipulations,  the  following  z- 
parameters  are  resolved: 


(39) 


2p 


1  1 


TTaa  (  1  -  0  1  2  ) 
o 


1  +  - 


tanh[<j)(0)-(f>(A)]  _ 


2p 


z  =  z 
12  2  1 


TTa[ao(ao+aX)]  2  (  1  -  0  1  2  ) 


nh  [(f)  (0)  -<J)  (A)] 


(3-54} 


“  2  p 


TTOt  (  a  +aA )  ( 1  -  0  1  2  ) 


1  - 


tanh£<t>(0)-(()(A)]  _ 


The  expressions  for  the  z-parameters  in  (3"54}  can  be 
altered  to  a  more  convenient  form  by  making  use  of  the 
following  relationships: 


2a 


i  _ 


TTCX 


0  where:  0  = 


a2A2  pcQA 

+  j  oo  - 


4a  2 

o 


it  a ' 


tanh (-x)  =  -  tanh  x 

cf>  (  A  )  -  4>  (0)  = 


aA 

10  1  1  n  (  1  +  — ) 


—  0  1  n  (  1  +  — ) 
OtA  a  n 


o  4pC 

i  -  e 1 2  =  -  j  oi  — 2. 

TTOt  2 


* 


' 


(4o) 


The  final  form  of  the  z-parameters  is  presented  as 
f o 1  lows : 


1  1 


j  LO  C  X 


0 


aX 


L  t  a  n  h  [—2.0  1  n  (  1  +— -)  ] 
aX  an 


2  a 


(3-55) 


Z  1  2  Z21 


.  aL  1  ra0  aX  _ 

J°JcoA(,+i7)  sinh[HX01n(1+— )] 


(3-56} 


1 


22 


j  (DC  X  (  l+~-) 
0  a0 


0 


aX 


tanh[^0  1n(1+^-)l  2ao 
aX  aQ 


(3-57) 


Let  : 


f  (ot)  =  ~  1  n  ( 1  +  2l) 
aX  a0 


(3-58} 


In  matrix  form:  V]  =  [z]  [l 


V 


V 


j0JcQX 


0 


'  ct'X 


1 


0 


tanh  f(a)0  2a 


(  1+— )  '2  s  i  nh  f  (a)  0 


1 


0 


(  1 +^_)  2  sinh  f(a)0  1  + 


aX 


0 


aX 


_tanh  f  ( a)  0  ^ao- 


(3-59) 


(41) 


In  the  limit,  as  a  approaches  zero: 


f  (a) 

0 


1 

[jay 


PcoX' 

"ao 


1 

2 


and  the  z-parameters  approach  those  derived  by  Castro 

and  Happ^^  for  a  distributed  R-C  network  with  no 

taper  (a  uniform  line). 

With  the  derivation  of  the  network's  z  - 

parameters,  these  can  be  used  in  the  subsequent  deriv- 

(33) 

ation  of  any  set  of  network  parameters  of  interest. 

For  example,  the  hybrid  or  h-parameters  may  be  derived 
on  the  following  basis: 


de  t  Z  z  1  2 

z22  z2  2 


'z2  1  1 
z2  2  z22 


* 


(42) 


3*3  Analysis  of  the  Distributed  R-C  Line 
Part  II  -  Network  Function  H(s) 

3-3-1  Introduction 


(30) 

In  1  965,  P r o t on o t a r I  os  and  Wing  presented  a 

significant  contribution  to  the  analysis  of  non-uniform 
distributed  R-C  lines.  Their  treatment  circumvents  the 
difficulties  associated  with  the  solution  of  the  govern¬ 
ing  differential  equation  and  provides  a  simplified 
method  of  deriving  the  voltage  transfer  network  func¬ 
tion  H  (  s )  .  The  network  function  H ( s )  is  extremely  use¬ 
ful  in  network  analysis  and  provides  information  re- 

04) 

garding  the  time  domain  behavior  of  the  line. 

In  their  approach,  the  distributed  network  is 
considered  as  a  chain  of  cascaded  two-ports  with  infini— 
tesmal  lumped  elements.  The  ABCD  parameters  are  derived 
for  a  single  infinitesmal  section  and  the  product  of  n 
elemental  transmission  matrices,  taken  in  the  limit  as 
n  -*•  00 ,  yields  the  transmission  parameters  of  the  whole 
line.  The  resulting  parameters  consist  of  uniformly 
convergent  series  which  are  entire  functions  of  s.  A 
summary  of  the  results  is  presented  in  Appendix  V, 

3.3.2  The  Network  Function  H ( s ) 

The  network  function  H(s)  can  be  derived  direct 
ly  from  the  transmission  parameter  A(s)  which  is 


defined  as  follows: 


■ 


' 


' 


(43) 


A  is  the  inverse  of  the  voltage  transfer  function 
from  port  1  to  port  2  with  port  2  open-circuited. 

Figure  I  I  I  -  3 

Transmission  Parameter  A 

For  considerations  regarding  an  open- c i rcu i ted  R-C 
line,  H  ( s )  can  be  simply  defined  as  follows: 

H  (  s  )  =  — — 

A  ( s ) 

3-3-3  Transmiss i on  Pa  rameter  A ( s ) 


From  Appendix  V 


A  ( s )  =  1  +  I  a  s 

n=  1  n 


(3-60} 


where: 


n  o 


X  x_  xo „ - 1 

/  c(x  )  /  2nr(x  ,)  f  c(x„ 

2n  o  Z  n  I  o 


2  n  -  2 


•  •  •  • 


(3-61 } 


/  3c(x_)  /  r{x})  d  x  ^  dx2  •  •  •  *  dx 


2  o 


2  n 


For  the  network  of  interest,  from  {2-1}  and  {2-2}: 


c  (x)  = 


c 

o 


{3-62} 


r  (x) 


P 

tt  (  a  +ax)  2 
o 


{3-63} 


Equations  {3-60,  { 3"  62 }  and  {3-63}  are  now  employed 

to  derive  the  coefficients  a^,  a^  and  a^. 


■  . 


(H) 


evaluation  of  a  ^ : 

From  equation  {  3  ~  6  1  } : 


X 


x, 


/  c(x0)  /  r ( x  )  dx  d x 

o  2  o  112 

X  x 

/  c (x)  /  r (y )  dy  dx 


(3-64} 


Substitution  of  equations  {  3  -  6  2  }  and  {  3  ~  6  3  3"  yields 
X  x  p 

a  =  f  c  /  — - - —  d  y  d  x 

1  o  o  o  tt  (  a  Q+ay )  z 


N  ow 


x  p  dy 


x 


o  ^  tt  (  aQ+ay ) 


it  a  (a  +ax) 
o  o 


Therefore: 


P  c Q  x  xdx 


f 


Pc0  r-  ^  1  _  ,  .  aX 


a  =  ' 

tt  a  o  a  +ax  7Ta 
o  o 


[  ~  ■  r  1  n  ( 1+~)  ] 


a 


(3-65) 


For  X  =  y : 

pc 


Vy)  =  ™ 


O  [  Ji 

L  a 

o 


1 


ay 


ln(1+-4  ] 
a  a_ 


{3-66} 


The  function  a^(y)  above  will  prove  useful  in  the  eval¬ 
uation  of  the  next  coefficient  a  £  •  Generally,  the 
integrations  to  follow  are  straight  forward  but  lengthy. 
For  the  sake  of  brevity  only  the  outlining  steps  are 


given. 


. 


eva luation  of  a ^ : 


(45) 


From  { 3  ~  6 1  }  : 

X  x  i  x  ^  Xo 

a  =  /  c  ( x I )  /  r ( x  )  /  c(x  )  /  r (x  ) dx  dx  dx  dx. 

2  o  4  o  3o  2o  1123  4 


A  change  in  the  bases  results  in: 


X  x  y  u 

a  =  /  c ( x )  /  r ( y )  /  c(u)  /  r(v)  dv  du  dy  dx 

2  0  0  0  0 

From  equation  {3-64}  and  {3_66}>  recognize: 


(3-67) 


y  u 

/  c  (  u )  /  r ( v )  dv  du  =  a  (y ) 

0  0  I 


Therefore: 


X  x 

a  =  /  c (x)  /  r (y ) a  (y )  dy  dx 

2  0  0 


{3-68} 


Designate: 


f (x)  =  /  r (y) a  (y)  dy 

0 


(3-69) 


Substituting  {3-69}  into  {3-68}  yields 


X 

a  =  /  c (x) f (x)  dx 

2  0 


{3-70} 


From  {3-63}  and  { 3 “ ^  6  }  : 


r (y ) a  (y  ) 


p2°°  (-  _1  v 
tt  2  ot  a0  (aQ+ay) 


1  ln(.1+^-) 


0+^)2 


■ 


(46) 


Completing  the  integration  in  {  3 “ 6  9 }  yields 


f(x)  = 


P2c, 


2  3 

it  a  a, 


/  axv 

1 n ( 1 + a^) 


,  / ,  aXN 

+  1  n  (  1  +— ) 

ao 


2ax 

a  +  ax 
o 


(3-71  ) 


Finally,  using  (3“7'1}  and  completing  the  integration 
in  { 3 “ 7  0  }  yields: 


2  2 
-  p  Co 

2  TT2a2 


,  o / ,  aXx  1  ,  „  aXN  .  , ,  ,  3X 

In  (l+~)  +—(3+ — )  1  n  (  1  + — )  - 


2a' 


a 


o 


aa 


The  substitution  of  y  for  X  in  equation  { 3 “ 7 2 } 
yields  the  function  a  £  ( y  )  which  will  be  useful  in 
evaluating  a  ^  • 


(3-72} 


(y)  = 


.  2  _  2 
P  cr> 


2  '  Tr2a2 


2a 


i  ay  1  ay  ay  3y 

1  n  2  (  1  +— — )  +  —(3+ — )  1  n  (  1  +~  )  -  — 


a 


3  oa 


e  v a  1 u  a  t i on  of  a ^  : 

From  {  3  _  6 1 } : 


X  x  x  g  x4 

f  c(*s)  f  6  r ( x 5 )  /  c(x4)  /  r(x3) 

°  b  0  0  v  0 

A  3  A  2 

/  c ( x 9 )  /  r (x  )  dxx dx2 •  • 

0  0 


•  •  • 


{3-73} 


dx, 


A  change  in  bases  results  in: 

X  x  y  u 

a_  =  /  c(x)  /  r ( y )  /  c(u)  /  r(v) 

3  0  0  v  0  g  0 

. f  c(g)  /  r(h)  dh  dg  dv  du  dy  dx  ( 3  “  7  ^ } 

0  ’  0 


0 


■ 


' 


(^7) 


From  { 3 “ 6  7  )  and  {  3  “  7  3  )  recognize: 

y  u  v  g 

/  c ( u )  /  r ( v )  /  c ( g )  /  r(h)dhdgdvdu  =  a  (y) 


The  re  f o  re 


X 


x 


a  =  /  c ( x )  /  r ( y ) a  ( y )  dydx 

3  o  o  2 


(3-75) 


Designate: 


g (x)  =  /  r  (y )  a  (y )  dy 

0  t- 


(3-76} 


Substituting  {  3 - 7  6  }  into  {  3  “  7  5  }  yields 


a  _  =  /  c(x)g(x)  dx 

i  0 


(3-77) 


From  {3-63}  and  {  3  ~  7  2 } : 


r  (y )  a  (y ) 


P3Cq 

3  2 

7T  a 


1  ln2Q^) 


/  ay  \ 

3  In  (1+74 


'SL 


2a2a2Q  (1+£V.)  2 

ao 

ay 

y  ln(1+l^)  3 


2  n  2 

o  (  1  +— ■) 
ao 


a“ a :  , .  ay, 2 


+  «a„  (1  +  “X)2 

3  o 


aa3  (i+^Vj 

o  '  +  aQ'  -1 


Completing  the  integration  in  { 3 “ 7  6  }  yields: 


g  ( x)  = 


3  2 

K  o 

tt3  a5  a 


1 

1 


,  2/<  ax\ 

1  n2  (  1  +  aT:) 


0+“*) 
a  _ 


-  3 


.  /  ,  otx. 

1 n ( 1 + ar) 


(1+^) 

a 

o 


.  ax,  .  6  ax 

3  n  <  +a  >  +  7"  777 

o  o  1  + — 


{3-78} 


> 


■ 


Finally,  using  { 3 “ 7  8 }  and  completing 


in  (3 


Note: 


(A8) 


-77)  yields: 


the  integration 


o  3  3 

P  c0 
tt3  a6 


[  -i  1  n  3  (  1  +— ) 


aX 


aX' 


(1+^)  1  n  2  (  1  +— — ) 


-  1  ln2(l+~)  -  b  (1+—)  ln(1+— ) 

a 


aX 

a 


aX 


6  1  n  ( 1  +~-~)  +  9  ~ 

a  a  _ 

o  o 


] 


{3-79} 


The  series  is  terminated  at  three  coefficients  at 
the  discretion  of  the  author.  It  is  felt  that 
sufficient  accuracy  is  attained  at  this  point 
and  any  extension  of  the  analysis  will  gain  very 


little  in  overall  accuracy. 


- 


' 

3.4  Extension 


(49) 

of  the  Analysis  to  Yield  the  Network  Function 
of  the  Complete  Equivalent 

3.4.1  General 

In  the  preceding  analysis  the  methods  of  P  rotono- 
tarios  and  Wing  were  employed  in  deriving  the  network 
function  of  the  distributed  R-C  network.  Since  the 
complete  equivalent  consists  of  the  R_C  network  shunted 
by  the  capacitance  Cj  an  extension  of  the  analysis  is 
necessary. 

In  the  analysis  to  follow,  modified  transmission 
parameters  are  derived  in  terms  of  the  transmission 
parameters  of  the  two  components  of  the  complete  equiv 

a  1  e  n  t  . 

3.4.2  The  Modified  Transmission  Parameters 


r  ~ 

r  (x)  1 

i 

1 

1 

- 1  \A 

1  • — 

1 

1 

t 

1 

/vAA - i 

-  1 

C  (  x )  | 

1 

- 

1 

1 

1 

C  i  1 

1 

1 

» - — - - 

1  1 

1  * 

■  .1 

1  1 

i 

network  1  network  2 

Figure  III-l 

The  Complete  Equivalent 

Let  the  transmission  parameter  matrix  of  network  1  be 


. 


1 


(50) 


and  let 
b  e  : 


the  transmission  parameter  matrix  of  network  2 


— 

— 

"  '  ' 

A2 

8  2 

1 

0 

C  2 

°2 

s  C  . 

1 

1 

The  transmission  parameter  matrix  of  the  system  will  be 
equal  to  the  matrix  product  of  the  transmission  para¬ 
meter  matrices  of  the  equivalent's  two  cascaded  com¬ 
ponents  : 


A  ' 

B  1 

A, 

B, 

1 

0 

A  1  +  s  C  .  B  1 

B, 

C  1 

D  1 

ci 

Di 

• 

s  C  . 

1 

1 

C 1 +s  C . D  ^ 

°i 

{3-80} 


Of  primary  interest  is  the  transmission  parameter  A' 
where: 


A 1  =  A 1  +  sCjB1 


{3-81} 

and  the  network  voltage  transfer  function  of  the  com¬ 
plete  equivalent  is  defined  as: 


H  1  { 5  )  =  — J—  <3-82} 

A  1  (s  ) 

From  equation  { 3 “ 8 1  }  it  is  obvious  that  an  evaluation 
of  A 1  ( s )  necessitates  the  derivation  of  the  trans¬ 
mission  parameter  B^s). 

From  Appendix  V: 

X  00 

B  ^  ( s )  =  J  r ( x) dx  +  I  b  ^  s n  (3-83) 


. 

* 


■ 


(51) 


where: 


b  =  /  r(x  -)  / 

n  o  2  n+  1  o 


x 


X2  n  +  1  ,  .  *2n  . 

c(x.  )  /  r ( x 

2  n  o 


3  x2 

/  c(x0)  /  r(xj  d  x  1  d  x  0 

0  ^  0  1  I  z 


2n  -  1 

•  d  x 


) 


2  n+  1 


(3-84} 


As  before  we  have: 


r  (x)  = 


•  p 


7T  ( aQ+ax) 


c  ( x)  =  c 


Let: 


X 

b  =  /  r  (  x)  dx 

o  o 


(3-85) 


evaluation  of  b ^ : 


X 


f 


pdx 


pX 


0  0  tt  ( a  +ax) 

o 


o 


{3-86} 


evaluation  of  b  ^ : 

F  rom  (3“ 84} : 


X 


/  r  ( x  _ )  /  3c(xJ  /  ^(xj  dx  dx  dx 


b  =  j  r  v  a  _  ;  j  j  j  1  v  ~  i 

1  0  3  0  2  o  1 


1  2  3 


A  change  in  bases  results  in 


X  x  y 

b  =  /  r ( x )  /  c ( y )  /  r ( u )  du  dy  dx 

1  o  o  o 


(3-87) 


From  {3-85}  recognize 


/  r(u)  du  =  b  (y) 


Therefore: 


X  x 

b  =  /  r ( x )  /  c ( y ) b  ( y )  dy  dx 

1  0  0  o 


0 


{3-88} 


■ 


(52) 


Designate: 


f  (x)  =  /  c (y) b  (y)  dy 

o  o 


Substituting  {3-89}  into  {3-88}  yields 


X 

b4  =  /  r(x)f(x)  dx 

1  0 


From  { 3" 83  }  and  { 3 - 8  6 } : 


(3-89) 


{3-90} 


c  (y )  b  (y ) 
o 


Pc, 


”0  (u^) 
a  o 


Completing  the  integration  in  { 3 “ 89 }  yields 


ftx)  =  -  md+as)  ] 


ax 


7T  a 


Therefore,  from  {3"83)  and  (3"92}: 
r (x) f (x)  = 


(3-91) 


(3-92} 


P2c( 


2  2 
7T  a 


ax 


/  ax. 

ln(1+— ) 


1  o 


a0(a0+ax) 


(  aQ+ax)  2- 


(3-93) 


Finally,  using  { 3  ~  9  3  ^  and  completing  the  integrati 
i n  (3-90}  yields: 


on 


P2c0 
^  tt  2  a 3  a 


,  1  n  (  1  +75—) 

1  n  (  1  +— )  +  - -  -  2 

ao  (1+sL) 

ao 


aX/aQ 

/  aX  \ 

1+T“) 


{3-9^} 


The  substitution  of  y  for  X  in  the  above  equation  gives 
the  function  b^(y)  which  will  be  useful  in  the  evalua¬ 


tion  of  b  2  '• 
P2co 

b,(y>  -  — — 

1  tt  a 3  a. 


/  aY  N 

ay  1  n  1  +  3  0 

1  n  (  1  + — )  +  - 

a°  (1+“£) 


ay  /a 


ay 
a  o 


2 


(3-95) 


n 


■ 


■ 


(53) 


evaluation  of  b2 : 

From  {  3  -  8  }  : 

X  x  5  x4  x3  x 2 

b  =  /  r(x5)  /  c(x4)  /  r(x3)  /  c(x2)  /  r(x.)*  •  • 

2  o  b  o  o  o  o  1 

.  dXj  dx2  dx3  dx4  dxg 


A  change  in  the  bases  results  in: 

X  x  y  u  v 

b  .  =  /  r  (  x)  /  c  ( y )  /  r(u)  /  c(v)  /  r  ( t )  d  t  d  v  •  •  •  d  x{  3  "9  6  } 

2  o  o  o  o  o 

From  {3-95}  recognize: 


y  u  v 

/  r(u)  /  c(v)  f  r(t)  dtdvdu  =  b  ( y ) 
0  0  1 


Therefore: 


X  x 

/  r ( x)  /  c (y ) b  1  (y )  dy  dx 
0 


(3-97) 


Des i gn  a  te  : 


g ( x)  =  /  c(y)b(y)  dy 

0  1 


(3-98} 


Substituting  (3“98}  into  { 3 “ 9  7 }  yields: 


X 

b  =  /  r ( x) g ( x)  dx 

2  0 


(3-99) 


From  {3-62}  and  (3-95): 


c (y ) b  (y )  = 

_  2 


^  22 
P  c  n 


7T  a0or 


ay 

, 

1  n  (  1  + — )  +  - 

a°  (1+HI) 


ay  /a 


-  2- 


(i+ssi) 

a  o 


(3-100} 


- 


' 


. 


(54) 


Completing  the  integration  in  { 3 “ 9  8 }  yields 


g  (x)  = 


_  2  2 
P  co 


2  4 

IT  a 


ax 


[  — 1  n2  (  1+— )  +  (1+— )  ln(1+— ) 


ax. 


+  2  1  n  ( 1+— )  -  3^21  ]  (3 


Therefore,  from  {3-63}  and  (3"101}: 


r (x) g (x)  = 


3  2 

P  co 

3  4  2 

tt  a  aQ 


1 

2 


r.  1n2(l+^)  ln(l+^) 


i  / i  ax\ 

1  n ( 1  +  0—) 


(1+— ) 

a 


axN  2 

o 


/  <  axN 
( 1+ — ) 

ao 


ax/a 


+  2- 


ax\  2 


(1+^) 

ao 


-  3- 


ax\  2 


( 1+«i) 

ao 


(3 


Finally,  using  (3-102}  and  completing  the  integrat 
in  (3-99)  yields: 


b2  “ 


3  2 

P  co 


3  5 

tt  a  a. 


aX 


1 


2 


i  2  /  i  Ot  X  \ 

in 

aX 

(1+T-) 


.  aX. 

In  1+— ) 
do 

,  ,  aX, 

(1+^ 


aX,  aX/ao 

3  ln(1+— )  +  - — 

(i+2i). 

ao 


-  101  } 


-102} 


i  on 


3 


(3-103} 


(55) 


3.4.3  Network  Function  H 1  ( s )  of  the  Complete _ Equi valent 

For  the  complete  equivalent,  by  previous  definition: 

H  1  (s)  =  — — 

A '  (s ) 

where: 

A1  (s)  =  A  +  s  C  j  B  1  {3-104} 

and: 

A  ^  =  1  +  a^s  +  a2s^  +  +  .  (3-105) 

B,  =  b  +  b  ,  s  +  b0s2  +  .  {3-106} 

1  o  1  2 

Substi  tut  ing  (3-105)  and  { 3 “ 1 0  6  }  into  { 3  -  1  0  4 }  yields: 

A'  (s)  =  1  +  { a 1  +  C  j  b  Q) s  +  (a2  +  C|b1)s2  +  (a^  +  C.b^s  +  . 

which  can  be  expressed  as: 

A 1  ( s )  =  1  +  a's  +  a's2  +  a  1  s  ^  +  .  (3-107) 

1  2  3 

where  the  modified  coefficients  are  defined  as  foil ow s 

a  1  =  a  „  +  C  .  b 

1  1  1  o 

a'  =  a  +  C.b  {3-108} 

2  2  '  1 

a '  =  a  +  C  b 

3  3  12 

In  general  form,  {3-107}  and  {3-108}  may  be  expressed 

A 1  ( s )  =  1  +  I  a’s"  {3-109} 

n=  1 

where: 

a 1  =  a  +  C . b  ,  {3-110} 

n  n  ' 


■ 

(56) 


and: 


H  1  (  s  )  = 


1 


1  +  l  a's 

n=  1  n 


(3-110 


The  expressions  for  coefficients  a  and  b  can  be  ex- 
r  n  n 

pressed  in  more  compact  form  by  using  the  following 


substi tutions: 


C  = 


pc 


o 


TTor 


n  =  1  + 


aX 


(3-1 12} 


6  = 


a 


c  a  ^ 
o  o 


Using  the  above  substitutions  yields 


a  i  =  C  (ri  -  1  -  1  nri )  (3-113} 

a 2  =  [i  1  n  2n+ (n  +  2)  1  nri"  3  (n- 1 )]  {  3  -  1  1  ^ } 

a ^  =  c3 [-iin3n  +  i(n-3) ln^n-^Cn  +  l) Inn  +  9(n-i)]  '(3-115} 

and: 


b  =  5  c 


n  - 1 

n 


6C2(=±ilnn-2^1) 


«C3  + 


{3-116} 

(3-117) 

{3-118} 


(57) 


3.5  Ne  two  rk  Function  of  the  Complete  System 

In  accordance  with  the  proposed  scheme  of  waveform 
restoration  (Figure  1-6),  the  complete  system  is  defined 
as  follows: 


Figure  I  I  I  -  5 
The  Complete  System 

The  network  function  for  the  complete  system  is  defined  as: 

V 

H"(s)  =  --  (3-1 19) 

Ve 

where:  V  =  restored  waveform 

o 

Vg  =  signal  applied  to  electrode 
k  =  open-loop  gain  of  amplifier 
With  perfect  restoration  H"(s)  =  1  for  all  frequencies. 

Note:  The  operational  amplifier  is  assumed  to  be  ideal  with 
infinite  input  impedance,  zero  output  impedance,  and 
open-loop  gain  k.  (see  Appendix  VI) 

The  development  to  follow  will  establish  the  frequency  de- 
pendance  of  H"(s)  and  is  based  on  a  knowledge  of  the  network 
function  of  the  complete  equivalent  H'(s). 


tm 


' 


■ 

- 


' 


(58) 


From  Figure  I  I  I  -  5  we  have: 

k  fVi  -  H '  (s) V  ]  =  kV,  -  kH '  (s) V 
L  1  o  '  o 


V  = 


This  yields 


Also: 


V 


V 


1  1  +  kH 1  (s) 


V 


V 


=  H  1  (  s  ) 


From  {3-119}: 

H"(s) 

From  {3-111}: 


V„  \J  ]  kH '  (s) 


V.  V 


V. 


1  +  kH '  (s) 


H 1  (s)  = 


1  +  a’s  +  a^s2  +  a^s3  + 


Combining  {  3 “ 1 2  0}  and  { 3 “  121 }  yields: 


where : 


H  "  (  s  )  = 


1 


k+ 1  1  +  q  5  +  q  5  2  +  0  s  3  4- 

2  3 


c  = 


1  k+1 


k+1 


c  = 


3  k+1 


{3-120} 


{3-121  } 


{3-122} 


{3-123} 


. 


(59) 


3  .  6  Network  Function  for  Woodbury's  Equalization 

(4) 

As  previously  stated  in  Chapter  I,  Woodbury  used 
restoration  by  equalization  based  on  the  simple  two-element 
model  of  the  electrode.  To  facilitate  a  comparison  of  the 
proposed  equalization  scheme  and  Woodbury's  first-order 
technique,  the  transfer  function  is  developed  for  the  com¬ 
plete  system  where  the  network  in  the  feedback  loop  of  the 
amplifier  is  the  simple  two-element  model.  This  system  is 
presented  In  Figure  MI-6. 


Figure  I  I  I  -  6 

Woodbury's  Equalization  System 

H'(s)  =  transfer  function  for  proposed  model 

R  =  total  electrode  resistance 

X 

=  /  r ( x) dx 

o 

C  =  total  shunting  capacitance 

X 

=  /  c (x) dx  +  C. 

o  1 

W(s)  =  transfer  function  of  V/oodbury's  model 


. 

' 


■ 


(60) 


The  network  function  for  Woodbury's  two-element  model  is  de¬ 
fined  as: 

V  9  1 

W(s)  =  —  =  —  ;  T=RC  {3-124} 

VG  UST 

For  the  feedback  amplifier: 

VQ  k  k( 1+sT) 

—  =  -  =  -  (3-125) 

V1  1  +  kW(s)  (  1 +k ) +s  T 

From  {3-111}  we  have: 

V1  1 

_  =  H  1  (  s )  =  -  {3-126} 

1  +  a' s  +  a ' s2  +  a ' s3  +• • 

C  12  3 


Using  {3-125}  and  {3-126}  the  overall  system  transfer 


t i on  be  comes : 

H1"  (s) 


VQ  VQ  V,  k( 1+sT) 
Ve  V,  Ve  ( 1  +  k) +s  T 


func- 


Th i s  may  be  expressed  as: 


k  1  +  sT 

H"  '  ( s )  =  - ; - - -  “  7* 

,  +  k  1  +  d ,  s  +  d2s2  +  d3s3  +  d^1*  + 


where: 


T 

+  - 

1  +  k 


+ 


a  IT 


1  +  k 


a  1  T 

d  =  a  1  +  ~  — 

3  3  1  +  k 


1  +  k 
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CHAPTER  IV 


NUMERI CAL  ANALYS I S 

b  .  1  1  n  t  r  o  d  u  c  t  i  on 

In  the  preceding  pages  of  Chapter  III  the  accepted 
equivalent  circuit  for  the  microelectrode  probe  in  its 
measuring  environment  was  analysed  mathematically  in  terms 
of  general  physical  constants. 

The  primary  objective  in  this  chapter  is  to  assume 
reasonable  numerical  values  for  the  describing  constants 
and  use  these  values  in  a  qualitative  evaluation  of  the 
complete  system's  stability.  The  relationships  presented 
in  section  3*5  will  form  the  basis  for  this  stability 
analysis. 

The  sections  to  follow  present  the  assumed  values  for 
the  describing  physical  constants  which  are  used  to  determine 
numerical  values  for  the  coefficients  characterizing  the  com¬ 
plete  system's  transfer  function.  The  numerical  analysis  is 
conducted  for  both  the  proposed  equalization  scheme  as  well 
as  for  Woodbury's  first-order  system.  A  computer  program 
used  for  Bode  plot  calculations  is  presented.  With  the  use 
of  the  Digital  Equipment  P  DP  8/S  computer.  Bode  plots  for  a 
number  of  test  cases  are  calculated.  These  results  are  used 
to  extract  information  regarding  the  frequency  and  time  do¬ 
main  ch a racte r i s t i cs  of  the  complete  system  proposed  for  the 
equalization  of  the  microelectrode. 
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4 . 2  Evaluation  of  Coefficients 
4.2.1  Numerical  Values  of  Constants 

Prior  to  the  evaluation  numerically  of  the  net¬ 
work  function  coefficients  assumed  values  are  required 
for  the  fundamental  constants  a,  X,  c  ,  a  and  p. 
taper  a 

A  value  of  1/60  is  assumed  in  accordance  with  Ling  and 
Gerard.^)  This  is  the  largest  tolerable  value  for 
satisfactory  results, 
line  length  X 

A  value  of  0.2  cm  is  assumed  in  accordance  with  Wood- 

.  (4) 

bury. 

t i p  radius  a 

A  value  of  1.67x10  ^  cm  is  assumed  based  on  a  tip 
diameter  of  0.5  microns  and  a  lumen  2/3  of  the  outer 
d i ame  te  r . 

capacitance  distribution  c  0 

A  value  of  10  pf/cm  is  assumed  based  on  the  results  of 
Appendix  VII  and  in  accordance  with  Nastuk  and 
Hodgkin.  ^ ) 

electrolyte  resistivity  p 

A  value  of  4  ohm-cm  Is  assumed  based  on  the  results  of 
Appendix  VIII. 
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4.2.2  Coefficients  for  the  Distributed  R-C  Line 

With  the  following  values  assumed  for  the 
amental  constants: 

p  =  4  ohm-cm 

a  =  1.67x10  3  cm 
o 

cQ  =  10  p  f /cm 
a  =  1/60 
X  =  0.2  cm 

substitution  in  to  { 3 -  1  12}  yields: 

£  =  4.58x10  8  Sec 
n  =  201 
6  =  1 01 **  fd“  1 

and  further  substitution  into  { 3 “ 1 1 3 1  (3"1l8} 

yields: 

-  6 

a  i  =  8.9x10  sec 
a2  =  0.98x10"12  sec2 
a  =  0.026xl0"18  sec3 

b  =  4 . 5x1 03  ohm 

b^  =  0.69  ohm2-fd 

b  2  =  3.8xl0~8  oh  m3 - f  d  2 

4.2.3  Coefficients  for  the  Complete  Equivalent 


fund- 


Using  the  numerical  values  established  in  sec- 
tion  4.2.2  and  the  relationships  of  { 3 “  1  0  8 }  yields  the 


. 
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a'  coefficients  for  the  complete  equivalent,  expressed 
n 

for  a  number  of  C.  values.  Numerical  values  for  these 
coefficients  are  presented  in  Table  I  V -  1 . 


Table  I  V -  1 


c 

(pf ) 

COE  FF 1 C 1  ENT 

2 

5 

1  0 

20 

a  |  ( x  1  0  ~  ^  sec) 

18 

31 

53.9 

98.9 

a^  (xIO  ^  2  s  e  c2  ) 

2  .  4 

4.4 

7.  88 

14.78 

a  ^  (xIO"  1  ^sec^) 

0.  1 

0.  22 

0 .  40 6 

0 .  786 

4.2.4  Coefficients  for  the  Complete  System 

Using  the  numerical  values  established  in  sec¬ 
tion  4.2.3  and  the  relationships  of  (3-123)  yields  the 

c  coefficients  for  the  complete  system.  Numerical 
n 

values  for  these  coefficients  are  presented  in 
Table  I  V- 2 . 

4.2.5  Coefficients  for  Woodbury's  System 

In  section  3*6  it  was  established  that  Wood 

bury 's  model  is  characterized  by: 

X 

R  =  /  r  (  x) dx 

0  X 

C  =  /  c ( x) dx  +  C . 

o 

Using  the  numerical  values  assumed  in  section  4.2.1 


we  have: 
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R  =  4.56  Mohms 

C  =  2  p  f  +  C  . 

r  1 

with:  T  =  RC 


we  h  a ve :  T  = 

18.24 

]isec 

for 

C  . 

1 

=  2 

pf 

T  = 

31.92 

ysec 

for 

c. 

1 

=  5 

pf 

T  = 

100.32 

ysec 

for 

c . 

1 

=  20 

pf 

Using  these  values,  the  numerical  values  established 

in  section  4.2.3  and  the  relationships  of  (3~127l 

yields  the  d  coefficients  for  Woodbury's  system. 

'  n 

Numerical  values  for  these  coefficients  are  presented 
in  Table  I  V  —  3  - 
4 . 3  Test  for  Stability 

With  numerical  values  established  for  the  coefficients 
of  the  complete  equivalent,  Routh's  stability  criterion 
(Appendix  I  X)  can  be  applied  to  the  characteristic  equation 
of  the  feedback  system: 


V. 


V 


a  ^  s 


+  as 
2 


+ 


+  (k+1) 


0 


' 
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Numerical  evaluation  of  the  Routhian  arrays  for  the  four 
test  cases  yields  the  following  information  regarding  system 
s  t  ab i  1 i  ty  : 

C.  (pf)  allowable  k  for  stability 

i 

2  -  1 <k<43  1 

5  -  1 <  k<  6  1  9 

10  -  1  <  k< 1 0  45 

20  -  1 <k< 1 859 

These  results  indicate  that  the  feedback  system  is  cond¬ 
itionally  stable  and  that  increasing  values  of  Cj  guarantee 
stability  for  wider  ranges  of  amplifier  open-loop  gain  k. 

An  examination  of  the  characteristic  equation  for 
Woodbury's  equalization  system  shows  that  this  system  is 
unconditionally  stable  for  all  positive  k. 

4 . 4  Program  for  Bode  Plot  Calculations 
4.4.1  Development 

For  the  general  case,  the  transfer  function 

takes  the  following  form:  00 

l  u  sn 

n  n 

,  U(s)  n=0 

G  (s)  =  -  =  - - 

Mis)  00 

v  v  s  y  m 

)  v  s 
m“n  m 
m=  u 

2 

U  +  UiS  +  UlS  +  •••••••••••••• 

s  )  =  - ■  - - -  — 

VQ  +  VjS  +  V2S2  +  . 

In  the  sinusoidal  steady-state  s  =  ja),  and  therefore: 


J 


I 


4" 
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( u  o~  u  2002  + . )  +  j(uia)-U2W3  +  *** . ) 

G(jaj)  =  - - - - - - - - - 

(  Vq“  v2^2+ . )  +  j  (  VjW- v3oj3+ . ) 

=  |  G  (  j  oo )  |  /  4>  ( oo ) 

where: 


(uq-u2oo2  + . )2  +  (u  1oo-u3oo3  +  *  •  • )  2 

|  G  (  j  oo)  |  = - - 

(  v0"v2l°2  + . )2  +  (  v  i  b)~  v3^3  +  *  *  * )  2 


and: 

4)  (oo) 


a  r  c  t  an 


u  .j  oo-  u  ^  oo  +  •  •  • 
u  -  u0oo2  +  •  •  •  • 

O  2 


V  1  00“  V300  3  +  •  •  •  • 

arctan  - 

v  - v  oo2  + . 

o  2 


For  the  transfer  function  of  interest: 

n  =  0 

m  =  3  . 

4 . 4 . 2  Computer  Program 

The  computer  program  to  follow  was  written  for 
the  case  whe  re : 

n  =  6 

m  =  7  • 

Provisions  are  made  for  the  arbitrary  selection  of  the 
desired  frequency  range  and  the  desired  sampling  den¬ 
sity  of  log-magnitude  calculations.  The  output  lists 
frequency,  absolute  magnitude,  magnitude  in  decibels 
and  phase  angle.  The  language  is  FORTRAN  IV  and  the 
machine  used  is  the  Digital  Equipment  P  DP  8/S  . 

C;  BODE  PLOTTER 

D  I  MENS  I  ON  A(6)  ,  B  (7)  lt  /N 

0;  F0RMAT(/,  "NUMERATOR  COEFFICIENTS  ARE:  ,/) 
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TYPE  0 
1  =  1 

1  ;  FORMAT (E) 

2  ;  ACCEPT  1 ,  A ( I ) 

1  =  1  +  1 

I F  (  I  -  7)  2,3,3 
3;  CONTINUE 

4;  FORMAT  (/ ,"DENOM  I  NATOR  COEFFICIENTS  ARE:11,/) 
TYPE  4 
1  =  1 

5;  ACCEPT  1  ,  B  (  I  ) 

1  =  1+1 

I F  ( I  -  8)  5,6,6 
6;  CONTINUE 

7;  FORMAT (/  /  'FREQUENCY (LOWER, FREQUENCY (HIGHER, 
DELTA" ,/) 

TYPE  7 

8;  FORMAT ( / , E , E , E , /) 

ACCEPT  8 , FI  ,  F2 , DF 
F=  F  1 

10;  W=F*2. 000*3. 14159 

RA= A ( 2 ) * ( W**  4 ) - A ( 4) * ( W**2 ) +A ( 6 ) 

AI=A( 1 )*(W**5)-A(3) *(W**3)+A(5)*W 
RB  =  -B  ( 1 ) *(W**6)+B (3) *(W**4) -B (5)*(W**2)+B (7) 
B  I  =  B  ( 2 ) * ( W**$ ) - B ( 4) * ( W** 3 ) +B ( 6 ) *W 
ABA=SQTF(RA*RA+AI *AI ) 

ABB=SQTF(RB*RB+BI *B I ) 

AN  G  A= ATN  F  (Al  /  RA) 

ANGA=ANGA*57. 295779 
I F  (  RA)  11,12,12 
11;  AN  G A=AN  G A+ 180.0 
12;  CONTINUE 

ANGB=ANGB*57 . 295779 

AN  G  B  =  ATN  F ( B I /RB) 

I F ( R  B )  13,14,14 

13;  AN  G  B  =  AN  G  B+ 1 80 . 0 
14;  CONTINUE 

GSV=ABA/ABB 
ANG.S  =  ANGA-ANGB 
DBS=20.0*LOGF(GSV)/2.303 
TYPE  1  6  ,  F  ,  GSV , DBS , ANGS 
16;  FORMAT ( /  ,  E  ,  E  ,  E  ,  E ) 

I F  (  F - F2 )  17,17,18 

17;  CONTINUE 

F=(10.0**(1 .0/DF))*F 
GO  TO  10 
18;  CONTINUE 
GO  TO  0 
END 
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4 . 5  Results  of  Bode  Plot  Calculations 

With  the  use  of  the  coefficients  of  sections  4.2.3, 
4.2.4  and  4.2.5,  the  computer  program  of  section  4.4.2,  and 
the  Digital  Equipment  PDP8/S  machine,  Bode  1 og - ma g n i t u de - 
phase  calculations  were  made  for  a  variety  of  test  cases. 
These  results  are  presented  graphically  in  Figures  IV-1  to 
IV- 3  in  the  form  of  magnitude  frequency  response  curves 
(fixed  C  j  ;  variable  k)  for  the  network  functions  outlined  in 
Table  I  V- 4 . 


Table  I  V- 4 


Tables  I  V- 5  and  I  V - 6  list  the  3db  (down)  cutoff  frequencies 
and  approximate  risetimes  for  each  of  the  selected  test  cases 
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Table  I  V- 5 

Unequal ized  Microelectrode:  H'(s) 


c . 

1 

f3db  (Hz) 

t  r *  (ysec) 

2 

r-'. 

oo 

40 

5 

5k 

70 

20 

1  .5k 

230 

* s e e  Appendix  X 

Table  I  V  -  6 

Equalized  Microelectrode 


C.  (pf) 

i 

open  -  1 oop  gain 

k 

PROPOSED 

WOODBURY  1 S 

H" 

( s  ) 

H"  1 

(s) 

f  3db 
(Hz) 

t  * 
r 

]i  s  e  c 

3  d  b 
(Hz) 

t  « 

r 

ysec 

2 

100 

1  .  5M 

.23 

780k 

.  45 

5 

1  00 

900k 

.  39 

500k 

.  70 

20 

100 

1  60k 

2  .  2 

1  60k 

2  .  2 

2 

200 

2 . 2M 

.  16 

1  .  1  7M 

.30 

5 

200 

1  .  7M 

.  2  1 

900k 

.39 

20 

200 

450k 

oo 

3  40k 

1  .03 

2 

500 

3  .  1  M 

.  1  1 

1  .  5  M 

.23 

5 

500 

2  .  5M 

.14 

1  .  3  M 

CNl 

20 

500 

1  .  4M 

.25 

700k 

.  50 
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k . 6  Conclusions 

Based  on  the  numerical  results  established  in  this 
chapter  the  following  may  be  concluded  with  respect  to: 
sy s  tern  stability  (section  4.3) 

For  the  proposed  equalization  scheme  the  system 
is  conditionally  stable  with  lower  and  upper  limits 
existing  for  the  amplifier  open-loop  gain  parameter  k 
(with  C.  fixed).  It  has  been  established  that  higher 
values  of  C.  guarantee  higher  allowable  values  of  k. 

For  Woodbury's  equalization  scheme  the  system  is 
unconditionally  stable  for  all  positive  k. 
low  frequency  (f<<f ^  ^ ^ )  error 

Since,  from  equations  (3-122}  and  { 3  “  1  2  7  }  - 

I  H"  ( s )  |  =  |H"'  (s)  |  =  -pk  for  f«f3db 

values  of  k  should  be  high  to  guarantee  minimum  magnitude 
errors. 

bandwidth  and  risetime 

On  the  basis  of  Bode  calculations  (Table  I  V -  6 )  , 
for  both  equalization  systems,  increasing  system  band¬ 
width  and  hence  decreasing  system  risetime  are  realized 
with  decreasing  values  of  C.  and/or  increasing  values 

of  k . 

A  comparative  examination  of  both  systems 


i n  d i cates  : 


I 


... 
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1.  A  very  significant  improvement  in  r  i  s e t i me 
relative  to  the  unequalized  m i c roe  1 ect rode . 

2.  A  definite  improvement  in  risetime  for  the 
proposed  scheme  relative  to  Woodbury's  scheme 
(at  fixed  k) .  Maximum  improvement  is  realized 
a  t  1 ow  C  j  values. 

i.e.  for  k=  1  0  0  and 


=  2  p  f  : 

32 % 

e  x  te  n  s 

i  on 

i  n 

BW 

=  5pf  : 

80% 

ex  ten  s 

i  on 

i  n 

BW 

=  2  0  p  f  : 

0% 

ex  ten  s 

i  on 

i  n 

BW 

Note:  The  numerical  analysis  in  this  chapter  has  been  based 
on  assumed  values  for  the  governing  physical  constants 
which  yield  a  total  electrode  resistance  in  the  order 
of  5  megohms.  Interpretation  of  the  numerical  results 
should  be  made  with  this  factor  in  mind.  The  con¬ 
clusions  are  based  only  on  the  relative  comparison  of 
the  network  transfer  functions  analysed  and  the  in¬ 
herent  trends  established. 
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CHAPTER  V 


CONCLUSIONS  AND  SUGGESTIONS 
FOR  FURTHER  RESEARCH 

5  .  1  Concl us i ons 

The  author  feels  there  are  two  main  contributions  made 
by  this  investigation.  Firstly,  the  analysis  of  the  dist¬ 
ributed  R-C  line  for  the  special  case  encountered  and,  sec¬ 
ondly,  the  equalization  proposal  for  glass  m i c roe  1 e c t rodes 
used  in  intracellular  action  potential  measurements. 

As  was  mentioned  previously,  analysis  of  distributed 
parameter  networks  has  received  the  attention  of  many  in¬ 
vestigators  in  recent  years.  The  solution,  in  closed  form, 
of  the  general  partial  differential  equations  describing  an 
a rb  i  t ra ry '  { r  (x)  and  c(x)  functions}  non-uniform  R-C  line  has 
not  been  found.  Detailed  analysis  has  been  isolated  to  a 
number  of  special  cases.  The  author  feels  that  this  investi¬ 
gation  provides  information  of  interest  to  these  researchers 
in  that  the  special  case  peculiar  to  this  thesis  is  unique. 

Of  particular  significance  is:  the  sinusoidal  steady-state 
solutions  for  line  voltage  and  current  (pages  31-32);  the 
general  solution  via  Legendre's  transformation  (page  31);  the 
network  z-parameters  as  a  basis  for  establishing  any  two-port 
parameter  set  (page^O);  and  the  derivation  of  the  network 
forward  transfer  function  H(s)  for  the  open  circuit  termina- 


. 


, 
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t i on  case  (page  **3)  and  H'(s)  for  the  capacitive  termination 
case  (page  50)  . 

In  this  investigation  a  proposal  utilizing  waveform 
restoration  by  equalization  has  been  made  to  correct  inherent 
distortion  contributed  by  the  low-pass  nature  of  the  glass 
m i c roe  1 e ct rode  used  in  the  measurement  of  cellular  action 
potentials.  Analysis  has  yielded:  a  proposed  equivalent  cir¬ 
cuit;  a  detailed  mathematical  description  of  the  microelect¬ 
rode  based  on  the  assumed  electrical  equivalent;  a  transfer 
function  description  of  the  proposed  equalization  scheme; 
and  a  numerical  evaluation  of  a  number  of  test  cases.  On 
the  basis  of  this  analysis  the  author  concludes  that  the 
proposed  restoration  scheme  is  indeed  feasible;  however,  de¬ 
pending  on  measurement  requirements,  the  degree  of  improve¬ 
ment  to  be  gained  over  Woodbury's  scheme  may  not  justify 
further  investigation.  Woodbury's  scheme  has  the  definite 
advantage  of  being  unconditionally  stable  (within  the  con 
straints  of  the  assumptions  made)  for  all  positive  k;  with 
increasing  k  being  consistent  with  increasing  bandwidth  and 
decreasing  low  frequency  error.  The  advantage  associated 
with  the  proposed  scheme  is  increased  system  bandwidth  (for 
low  Cj)  due  to  more  exact  equalization  and  the  possibility 
of  utilizing  the  positive  boost  which  occurs  as  the  system 


approaches  instability. 
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5  .  2  Suggestions  for  Further  Research 

It  is  considered  beyond  the  realm  of  this  investiga¬ 
tion  to  make  specific  quantitative  design  proposals  which 
would  result  in  an  immediate  physical  realization  of  the 
restoration  circuitry.  However,  on  the  basis  of  the  in¬ 
vestigation  as  it  stands,  and  assuming  that  sufficient  motiv¬ 
ation  exists  for  improvements  to  the  degree  implied  by  the 
results  of  this  study,  further  research  could  be  directed  in 
the  following  areas: 

1.  A  numerical  analysis  for  a  wide  range  of  realizable 
electrodes  in  terms  of  assumed  physical  constants  to 
determine  whether  the  success  of  the  proposed  scheme 
is  electrode  dependent. 

2.  An  investigation  of  the  noise  characteristics  of 
the  proposed  equalization  scheme. 

3.  A  study  of  the  effect  of  non-ideal  operational 
amplifier  ch a rac te r i s t i cs  on  the  response  of  the  prop¬ 
osed  equalization  system. 

h.  A  study  of  the  possibility  of  economic  fabrication 
of  the  required  distributed  parameter  network  with 
provision  for  external  control  of  effective  line 

length. 
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APPENDIX  I 

DERIVATION  OF  FUNCTION  r(x)  FOR  THE  MICROELECTRODE 


Figure  A  I  -  1 


x  x+Ax 

Figure  A I  - 1  represents  a  differential  element  of  the 
electrolytic  column  with  a(x)  representing  the  radius  of  the 
column  at  position  x  along  the  column.  For  this  element: 


a ( x)  =  a  +ax 
o 

A  (  x)  =  7T  (  aQ+ax)  2 
a  ( x+Ax)  =  aQ+a(x+Ax) 

A  (x+Ax)  =  tt{  aQ+a  (  x+Ax)  }  2 

Consider  the  resistance  R  of  this  element: 

.  .  Ax  . 

AR  =  p  - 

average  area 


where:  average  area  =  i( A ( x) +A ( x+Ax) } 

=  A  (x)  +2-7ra(  2aQAx+2axAx  +  aAx2) 


Therefore: 

r  (x) 


=  1  i  m 

Ax->0 


AR 

Ax 


P 

MAO 


p 


dR 
d  x 


7T  ( a Q+ax)  2 
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APPEND  I  X  I  I 

DERIVATION  OF  THE  PARTIAL  DIFFERENTIAL  EQUATION  WHICH 
CHARACTERIZES  A  GENERAL  DISTRIBUTED  R-C  NETWORK^12) 


I  (x  ,  t ) 


A 

-9- 


V (x  ,  t) 


r  ( x)  Ax 

A/WA- 


i  ( x  +  Ax  ,  t ) 
c ( x) Ax 


r ( x+Ax) Ax 

— WA/VA 


A  i  (x  ,  t) 


v ( x+Ax  ,  t) 


0 


c ( x+Ax) Ax 


x+A  x 


Figure  All-2 
A  Differential  Element 

The  foil ow  ing  development  is  general  and  applicable  to  any 
arbitrarily  tapered  R-C  line.  Attention  is  restricted  to 
one  dimensional  current  flow.  The  introduction  of  tapers 
results  in  resistance  per  unit  length  r(x)  and  capacitance 
per  unit  length  c(x)  where  x  is  the  space  parameter  in  the 
direction  of  current  flow.  Figure  All-1  illustrates  the 


. 
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accepted  symbolic  representation  of  a  tapered  distributed 
R-C  network.  Figure  All-2  represents  a  differential  element 
of  this  network. 

Consider  the  loop  equation  of  ABO: 

v(x,t)  -  i  ( x+Ax  ,  t )  r ( x) Ax  -  v(x+Ax,t)  =  0 
Therefore: 

v  (x+Ax',  t)  -  v  (  x  ,  t ) 


Ax 


=  -  i ( x+Ax , t ) r ( x) 


Consider  the  limit  as  Ax+0;  the  equation  becomes: 


3v  (x  ’  t)-  =  -  i  ( x  ,  t)  r  (x) 


3  x 

Consider  the  node  equation  at  A: 

i  ( x , t )  =  A i  ( x  ,  t )  +  i  (x+Ax  ,  t) 


Also: 


A  i  ( x  ,  t )  =  c(x)Ax 


3  v  (  x  ,  t ) 

3 1 


{1} 


Therefore: 


(x,t)  =  c  ( x )  A  x  }v(x»t)-  +  i(x+Ax,t) 

3t 


This  becomes: 


\  (x+Ax, t)  -  ?  ( x , t) 

Ax 


-  c 


(x) 


3 v  ( x  ,  t) 

'  3tT 


Consider  the  limit  as  Ax->0  ;  the  equation  becomes 


3 i  (x  ,  t ) 
3x 


=  -  c 


(x)  3 v  ( x 

3 1 


{2} 
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'With  v  and  i  functions  of  x  and  t,  equations  {1}  and  {2} 
can  be  written: 


9  v 
9x 


=  -  r ( x) 


and 


From  {  1  }  : 


= 

9x 

1  9  v 


/  \  9  v 

c  ( x)  — 

9 1 


r ( x)  9x 


Substituting  into  {2}  yields 


8 
9  x 


1  9  v  1 


r  ( x)  9 x 


=  -  c 


(x) 


9  v 


9 1 


Expanding  the  partial  differentiation  yields: 


1  9  z  v 


1  d  r  (  x)  9  v 


r  ( x )  9 x 2  {  r  (  x)  }  2  dx 


9  x 


( x) 


9  v 

9  t 


Finally  the  equation  becomes: 


9 2  v  1  d  r ( x)  9  v  9v 

-  - - -  r  (  x)  c  ( x)  — 

3x2  r ( x )  dx  9x  9t 


A  similar  development  for  current  yields: 


(3) 


9 2  i  1  d  c ( x)  9  i  . 

-  - - =  r  (  x)  c  (x) 

9x2  c(x)  dx  9x 


{4} 


Both  equations  (3>  and  {4}  govern  the  behavior  of  a  dist¬ 
ributed  R-C  line  in  which  R  and  C  are  functions  of  x  only. 


. 
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APPEND  I  XIII 


TRANSFORMATION  OF  EQUATION  {3-1*0 


Equation  {3“1*0  has  been  expressed  as 


(a^+ax)  2  - — 2a(a^+ax)  ~  +  — —  k2X  =  0 

d  x 


dx 


o 


{1} 


TT 


This  is  an  ordinary  linear  second-order  differential  equa¬ 
tion  with  variable  coefficients  and  is  recognized  as  the 

(  32 ) 

second-order  form  of  Legendre's  Linear  Equation.  The 

equation  can  be  reduced  to  a  linear  equation  with  constant 
coefficients  by  means  of  the  following  transformation: 

e  z  =  ( aQ+ax)  { 2 } 

The  evolution  of  the  transformation  is  as  follows: 

eZ  dz  =  a  dx 


a 


d  z 


dx  aQ+ax 


dX 
d  x 


d  2  z 


dx 

d  X  d  z 
d  z  dx 


a' 


( aQ+ax) 


a 


dX 


a0+ax  dz 


The  re  f o  re : 


,  •  dX  dX 

( a  +ax)  —  =  a  — 

dx  dz 


{  3  ) 


F u  r  th e  r  : 


d  2  X 


dx' 


dx 


d  dXdz  dz  ddX  dzdX 

-  [ - ]  = - -  [  —  ]  + - 

dz  dx  dx  dx  dz  dx2  dz 


. 


a 


dx  a  +ax  dz 
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d  2  X 


dx' 


•a' 


c 


d2X  dX 


(aft+ax)  dz  dz 


Therefore: 


(a  +ax) 
o 


d2X 


d  x  ‘ 


=  a' 


[ 


d  2  X  dX 


dz2  dz 


] 


{4} 


Substituting  equations  (3)  and  {4}  into  equation  {1}  yields 
the  transformed  equation  (3)  where: 


X  =  X  (  z) 


X  (  x)  X  (  z) 


d2X  dX  k  c  p 

-  +  —  +  - X  =  0 

dz2  dz  Trot 2 


(5) 
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APPENDIX  IV 


CONDITIONS  ESTABLISHING  det[A]=0 


Consider  the  following: 


[B]  =  [A] 


-  1 


de  t  [A] 


1  j-  a  22 

_a2  1 


-  a 


12 


1 


1 1 


where:  det[A]  -  a]]a22  ~  a  1  2 a 2  1 

From  previously  developed  expressions: 


de  t  [a]  =  -  (^)2{ao(ao  +  aA)}2(1-0,2)sinh{^(O)-4)(A)} 

From  the  above  expression: 

1.  de  t  [  A] =0  for  (  1  -  0 '  2 ) =  0  or  0,2=1 

2.  det[A]=0  for  s  i  n  h  {  cf)  (  0 )  -  4>  (  A )  }  =  0  or  cj)  (  0  )  -  (j)  (  X  )  =  0 
Expanding  on  these  statements: 


statement  1: 


0 


I  — 


i  ,  •  4pCo 

1  +  J  00 - j 

ma 


0  c 

0  '  2  =  1  f  o  r  j  co — — %  =0  or  0)  =  0 


not 


statement  2 


conclusion: 


(<j)  (0)  -  <J)  (X)  }  =  £6  1  1  n  ( 


aQ+aX 


-)  =  0 


on 


ly  possible  for 


£_  =1  or  aX  =  0 


This  condition  is  not  applicable 


on 


which  will  yield  de  t  [a]  -  0  is 


co 


The  only  conditi 


0 


UV- 


. 
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APPENDIX  V 

ANALYSIS  OF  P ROTON OT AR I  OS  AND  W I N  G  ^ 3  °  ^ 

A  5  -  1  Transmission  Parameters  (ABCD) 

OO 

A  =  1  +  I  a  s n 
n=  1  n 

X  00 

B  =  /  r(x)dx  +  £  h  s 

0  n=  1  " 

X  00  n 

C  =  s  /  c(x) dx  +  s  £  cs 
o  i  n 

u  n  =  1 

00 

D  =  1  +  l  d  sn 
n=  1  " 

where: 

■  ,^c(x2n)„/X2nr(x2n-1) . Q /"  3  o  (  x  2 )  q /"  2  r  ( x  , '  d  x ,  •  -  d  x  2  n 

bn  -  /rU2n  +  ))/Zn+]cU2n) . (/3c(x2)(/X2r(x1)dxrdx2n+1 

cn  -  /^2n^/2n+'r{*ln] . o/X3r(x2).;X2c(X1)d><1"dX2n^ 

dn  -  o/r(x2„)o/X2nc(x2n-,) . /X3r(x2)0/X2c(x,)dx,--dx2n 

n  =  1,  2,  3,  . 

r(x)  =  resistance  per  unit  length  variation  function 
c  ( x )  =  capacitance  per  unit  length  variation  function 

=  length  of  line 


X 


(9*0 


The  four  series  are  all  uniformly  convergent  and  the  ABCD 
parameters  are  entire  functions  of  s. 

A  5  ~  2  The  Network  Function  H(s) 

From  the  definition  of  the  transmission  parameter  A; 


A  ( s )  1  +  alS  +  a2s2  +  a^3  +  . 

The  network  function  is  stable  with  poles  in  the  left  hand 
plane  only.  Also,  the  poles  of  the  network  function  of  an 
arbitrary  n  on- un i form  R  -  C  line  terminated  in  an  R~C  imped- 
ance  (including  open  and  short  circuits)  are  all  negative 
real. 

A  5 -  3  Time  Response  -  Delay  and  Rise  Time 

From  the  theory  developed  by  Elmore^^  the  delay  and  rise 
time  of  the  step  response  can  be  determined  from  the  coef¬ 
ficients  of  s  in  the  network  function  H(s).  For  a  tapered 

R  -  C  line: 

td  ■  ai 

TR  “  {2TT(a^-2a2) 


. 

* 
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APPENDIX  VI 


THE  OPERATIONAL  AMPLIFIER 


(3  7) 


The  operational  amplifier  Is  a  high  gain  d-c  amplifier 
ch a rac te r i zed  by:  -  near  infinite  gain 

-  near  infinite  input  impedance 

v  -  near  zero  output  impedance 

-  wide  frequency  bandwidth 

Consider  the  following  general  voltage  operational  amplifier 


where: 


s 


and  Z 
is  the 
is  the 
i  s  the 


are  short-circuit  t 
input  impedance 
output  impedance 
load  impedance 

( s ) 

- — -  is  the  transfer 

( s  ) 


amplifier  with 


rans  fe  r  i mpedan  ces 


function  of  the 
no  load 


. 


' 
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Circuit  analysis  yields  the  following  relationship: 


i  .  1  (i  +  1  +  '  > 

Zf  G  Z1  Zf  zi 


£°(i  +  l  +  i)i 

G  Z  If  Z.  o 


.  £l 


-  I 


For  the  conditions:  G ( s )  =  -  k 


Z  = 


I  =  0 
9 


the  relationship  becomes: 


Lif+  k  (i.+  i+  ; 1 


1 


With  k  very  large  and  Z  ^  ,  Z^.,  Z.  the  same  order  of  magnitude 

the  relationship  reduces  to: 


Eo(s)  Zf(s) 

_ _  a/  . 


E  1  (  s  ) 


Z  1  (  s  ) 


The  assumptions  which  must  hold  for  good  operational 
amplifier  design  are:  -  G(s)  =  -  k  and  k  large 

I  =  0 

g 

Z  very  small 
o 

Z.  large  and  in  the  order  of 
Z  ^  and  Z  ^ 
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APPENDI XVII 

NUME  R I  CAL  EVALUAT ION  OF  c 

- —  Q 

A  discussion  of  the  capacitance  distribution  in  section  2-3 

established  that  c(x)  was  constant  and: 

c  (x)  =  c  =  2tte/  1  n  (  r/a) 

o 

Numerically,  assuming  that  the  glass  used  is  Corning  7060 

Pyrex  (e  =  43.8xl0"12  fd/m)  and  that  the  ratio  r/a  is  1.5, 

we  find  the  capacitance  cq  to  be  0.68  p f /mm .  This  is  in 

(2) 

agreement  with  calculations  made  by  Nastuk  and  Hodgkin. 

The  value  assumed  for  this  investigation  will  be  1  pf/mm 
or  10  p  f /cm . 


■ 


■ 
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APPEND  I  XVIII 

NUMERICAL  EVALUATION  OF  p  (  3  M  K  C  1  ) 

Of  specific  interest  here  is  the  resistivity  of  the  elect¬ 
rolyte  ( 3 M KC 1 )  which  forms  the  conductive  path  within  the 
microelectrode.  The  following  information  was  extracted 
from  the  International  Critical  Tables  (Volume  Vl,page  229) 

electrolyte :  3MKCI 

(223.68  grams  dissolved  in  1  litre 

of  water) 

concentration :  C=3 0 0 0  milli-mole s /litre 

temperature:  18  degrees  C. 

106  K/C  =  88.1 

K  =  specific  conductance  (ohm  ^  - cm  ) 

From  the  above  information: 

88.1  x  3000  -  1  - 1 

K  =  -  =  0.264  ohm  -cm 

,06 

1 

Therefore:  P  =  =  3*78  oh m- cm  @  18  C. 

K 

In  this  investigation  the  resistivity  is  assumed  to  be 


4  ohm-cm. 


- 

. 
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APPEND  I  X  I  X 


APPLICATION  OF  ROUTH 1 S  STAB  I  L I TY  CRITERION 


(38) 


The  application  of  a  Routhian  array  based  on  the  characteris¬ 
tic  equation  of  the  closed  loop  transfer  function  of  a  feed¬ 
back  system  yields  information  regarding  system  stability. 


The  characteristic  equation  is  recognized  as 

1  +  k  H  1  (s)  =  0 


From  {3-11}: 


H  '  (  s  )  = 


1 


1  +  a's  +  a  '  s  z  +  a's 

1  2  3 


3 


This  yields  the  characteristic  equation 


3  2 
a's  +  a's 
3  2 


+  a  j  s  +  ( k+  1  )  =  0 


' 

. 
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The  Routhian  array  becomes 


a2a  1  "  a3  (k+1  ) 


( k+  1 ) 


0 


( k+  1  ) 


0 


0 


The  Routh  stability  criterion  can  be  stated  as  follows: 

"A  system  is  stable  if  the  Routhian  array  of  its 
characteristic  equation  contains  no  negative 
elements  in  its  first  column." 

With  k  a  definite  part  of  the  Routhian  array,  the  Routh 

dicate  the  range  of  k  which  will 


criterion  can  be  used  to  in 
give  s  tab  1 e  response . 


* 

- 

■ 


. 
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APPENDIX  X 


STEP  RESPONSE  OF  AN  AMPLIFIER 


(39) 


Under  the  conditions  that  the  input  to  an  amplifier  can  be 
represented  by  a  single  time  constant  low-pass  network,  it 
is  possible  to  relate  the  system's  rise  time  to  the  upper 
cutoff  frequency  in  the  following  manner: 


T=- 


RC. 


1 _ 

2. -nr  T 


2.2  0.35 

t  =  2  .  2  T  = - T -  =  - 

r  2  71  f  3  d  b  f  3  d  b 

With  systems  which  are  not  single  time  constant  in  nature  the 
above  relationship  can  be  used  as  a  first  approximation  of 


system  response. 


: 


> 


